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On the Size-Frequency Distribution of Penetrating Showers 


By W. HEITLER anp L. JANOSSY 
Dublin Institute for Advanced Studies 


MS. received 30th Fune 1949 


ABSTRACT. We have calculated the relative number of showers containing 2 mesons 
when a fast primary nucleon hits a compound nucleus on the assumption that in a single 
nucleon-nucleon collision only one meson is produced (pure plural production). The 
cross section for meson production is assumed to have the form (¢/E) de/E, where E is the 
energy of the primary and « the energy lost by it, but the absolute value of ® and its depen- 
dence on </E are left open. ‘The fluctuations of the shower size are duly taken into account. 
For comparison with experiments it is assumed that about half the number of fast shower 
particles are fast protons as suggested by the positive excess of penetrating shower particles. 
The theoretical results depend on two constants, the total cross section and some weighted 
average over the energy distribution ®(</E). Choosing these suitably, the results are in 
good agreement with the relevant data within the experimental statistical errors, and the 
agreement extends also to the “‘ large showers ”’ with 20 or so thin tracks. The values of 
the constants thus determined are physically reasonable. It is concluded that the size— 
frequency distribution of penetrating showers can be fully understood on the basis of pure 
plural production, but some—small or large—contribution from multiple processes cannot 
be excluded. 


§1. INTRODUCTION 


Nrecent papers by Brown, Camerini, Fowler, Heitler, King and Powell (1949) 
|e relative frequency of stars with a given number of ‘thin tracks” is 
measured. ‘There can be little doubt that these events in the photographic 
emulsion are identical with the penetrating showers previously found and 
observed in counter and cloud chamber experiments by Janossy, Rochester and 
others. Ifthisis true, then it is also highly probable that we have to deal here with 
the process of meson-production by fast nucleons (and, indeed, roughly half of 
the stars are initiated by charged, half by neutral particles) and that at least a 
considerable fraction of the thin tracks are mesons (z-mesons, presumably, 
perhaps also 7-mesons). Some of the tracks are also likely to be fast protons. 
For the interpretation of these events, two widely different hypotheses have 
been put forward: (i) The occurrence of groups of particles is explained in a 
natural way as due to a multiple collison of the primary nucleon with several of the 
closely packed nucleons in a compound nucleus, O, C, N or Ag, Br in the photo- 
graphic plate. We call this process, following J. G. Wilson, plural production. 
(ii) The groups are assigned to a single elementary process in which the fast 
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nucleon would produce several mesons in one single collision with a single nuclear 
nucleon (genuine multiple process). It is assumed here that the emission of 
several mesons in one act is more probable than that of a single meson. The 
process (ii) would be a type of process for which so far no analogue is known in 
physics*, in so far as the emission of the minimum number of particles or quanta 
(compatible with conservation laws) is always the most probable event. A 
decision between the two interpretations therefore involves a rather fundamental 
issue and would be of great importance for the future development of meson 
theory. In this paper we examine the question whether the data available to 
date can be understood by the assumption (i) of pure plural production. It will 
be seen that this is the case quantitatively within the present accuracy of the experi- 
mental data, with exception of an irregularity that—possibly—occurs in the 
experimental distribution curve (cf. §6). But this irregularity is not firmly 
established. 

The agreement also applies to the ‘‘large showers.” observed so far (up to 
22 particles and, perhaps, even to the one case of 31 particles, cf. §6), and there is 
no necessity to assume a different type of process for these showers. 

It will furthermore be seen that the assumptions one has to make about the 
cross section for meson production etc. to reach this agreement are physically 
reasonable and they are in harmony with the known facts about the absorption 
coefficient of fast nucleonst. ‘They are further compatible with the values given 
in the theory of meson production put forward by Hamilton, Heitler and Peng in 
its later version}, within the rather wide limits set by this theory. 

Whilst the facts can be explained by the assumption of pure plural production, 
the issue is not yet decided. It might be possible to explain the facts also by 
process (ii) (see, e.g., Heisenberg 1949, also Dallaporta and Clementel 1948), 
making suitable assumptions about the multiplicities and their energy dependence, 
but a complete theory that would predict them as more probable than the 
single emission is not available. Nor can it be stated—even if the explanation 
of mainly plural processes is accepted—that all of the showers in question are 
plural. It might well be that a fraction of the showers or part of the multiplicity 
in a given shower is due to genuine multiple processes. 

We shall work out the pluralities of the process in a more phenomenological 
way without assuming any detailed theory of meson production. We shall only 
make some very general assumptions about the energy dependence of the cross 
section etc. It will be seen that then two constants occur, the total cross section 
and some weighted average over the energy distribution of the emitted mesons. 
We shall try to fix these constants, as far as is possible, by comparison with the 
facts. 

For a first orientation we shall neglect in this paper any secondary production 
of mesons by fast recoil nucleons within the same nucleus. This effect may not 
be negligible when the primary is very fast and tends to make an already big shower 
still bigger, especially in heavy nuclei. This will have to be taken into account 
in the discussion of the available material (§ 6). 


> 


* The emission of several extremely soft photons in any electromagnetic process (which is 
sometimes quoted as an analogue) hardly bears any resemblance to the high energy multiple 
processes in question. ‘The former is a purely classical phenomenon, translated into the language 
of quantum theory. 

t Cf. Heitler and Janossy 1949, quoted in the following as HJ. 

t For references see Heitler 1949, quoted in the following as HHP. 
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§2. THEORY OF PLURAL PRODUCTION 


We wish to calculate the probability for mesons to be emitted in a collision 
of a fast nucleon with a given nucleus. Let ®(e, Z)de be the cross section for 


7 


emission of a meson with an energy loss « suffered by the nucleon with initial 
energy FE in a collision with a single nucleon. « includes the recoil energy. 


Instead of « we use E’=E-—c. Ina passage through a nucleus in general several 


-acts of meson emission and discontinuous energy losses will occur. We make the 
following basic assumptions : 


(i) "That O(E’, E) dE’ =® (=) as | “MT (1) 

SORIA 
shall bea function of E’/E only. ‘This assumption has a high degree of probability 
for being at least approximately correct. In the first place the analogous law 
for Bremsstrahlung has the same form, and in the second place the (very crude) 
theory of HHP leads to the same result. (1) holds, of course, only for sufficiently 
high energies E, say above a certain critical energy E,. 

(ii) ‘The individual acts of emission shall be Seaieteitis independent events. 
‘This assumption is somewhat doubtful, on account of the large disturbance 
caused in the nucleus even after the first emission. However, one may argue that 
‘the recoil energies produced in each emission act are in most cases small compared 
with the energy of the primary nucleon. The energy communicated through the 
recoil to the second nucleon which suffers the second hit from the primary will 
be even smaller, and since we may treat the nuclear nucleons anyhow as free it 
will hardly be very wrong to consider them all at rest before they are hit by the 
primary. 

It is convenient to write 


E’\ dE’ E’\ dE’ Be LEN dE: 
o(5)> =0.0(5) >. {o(s)F-! oh iis) 9 tale (2) 


(w differs from the quantity w used in HJ by the factor Nd,®). © is the total 
cross section and is, by (1), a constant. Let N be the density of nuclear matter. 
‘The probability for a total loss of energy from E to E’=E,,_, in n—1 individual acts 
of emission, each taking place in the intervals of distance between x, and 
Net Ova. X,-4 abd X, 4 +aX, 45 18 


pte fo fo BB o(B)E 


E dE 
kot —— ) 1 (ON) exp (—ONx,,)dx,.:..d%p_y) -+0-- (3) 
Esa n—2 
es Ce Ob hs ren 0 al ae OUI ed ORs Ae Raa bs ara (Gr) 
'The integral is to be extended over all values £,, E,....E,~2 (Ey-1 18 given) 
satisfying (3’). It is furthermore to be extended over all ie OPM Na ak rei ty 
satisfying 
es PNM pM SAV | iene se ols (37) 
"The factor exp(—@®Nx,,) is the probability that within the total distance w,, no 
emission act takes place, except in the infinitesimal intervals dx,....dx,_, whose 


total length is, of course, zero compared with x,. ‘The reason for denotiny the 


itotal length by x,, will become clear below. 


We introduce V_ =log (Em_1/Em), Y =log (£/E’). 
46-2 
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The integration over x,....%,_1 yields «,,”"1/(m—1)! because the integration can. 
be performed for each x,, from 0 to x, irrespective of order, and keeping the order 
(3”) introduces merely the factor 1/(n—1)!. Then (3) becomes 


Pn1dE' = [eGror) oe ODn1)\OVy Va sae OV; 29 


I<) page ae 4) 
dE - (ONx,)" ( 
Ei exp (—®Nx,) arrier ye 


x 


yD Ve T= nel = te e (4’) 


and the normalization for w becomes 


J, 0m) dyn =1. 


The length of path in nuclear matter is, of course, limited, and clearly two 
cases are possible. (1) When emission acts take place the primary nucleou 
may have lost so much energy that it is no longer capable of further meson emission. 
This is the case if, in the last emission act, the mth say, the nucleon loses energy 
from a value E’>£E, to some value below £,. (11) The primary nucleon may 
commit 2 emission acts and leave the nucleus still retaining an energy >E£.. 
We call the probabilities for these two types of events p® and p® respectively. 

We consider first the case where one further collision takes place, the mth, 
in which E drops from £” to a value E” below the critical energy F,. Below E,, 
w will steadily decrease and will eventually vanish at the meson threshold. We 
have no exact information as to how rapid this decrease is, but there are reasons, 
both experimental and theoretical, for believing that it is fairly rapid. We can 
idealize the situation by assuming that below EF, the probability for further meson 
production can be neglected, but that for the jump from £’ to E” the cross section 
has still its full value w®. Evidently this idealization is crude, but since it is 
made only for the last step the total error cannot be very large, except for small 
showers with only one or two particles which we shall exclude from the comparison 
with the experiments. The probability for the last step is 


Bek (Ee Nan 
on | w (F) Ft 
Instead of carrying out the integration over the y,,’s in (4) we first average over 
the primary spectrum of the incident nucleons. In HJ it was shown that this 
remains a power spectrum for any depth in the atmosphere if the spectrum falling 
on the top of the atmosphere was a power spectrum, and this is very probable. 
We therefore assume a primary spectrum of the form 


F(E\E =y Ey 


Fa Doe SS SR ee (5) 
(5) is normalized so that the total number of nucleons above the critical energy F, 
isunity. Primaries with E<£, do not contribute according to our assumptions. 
Actually (5) is only correct if #, is larger than the latitude cut-off energy (for 
Europe 3 x 10° ev.). We have reasons to think that FE, is of about the same order 
of magnitude (again the theory of HHP suggests this). Even for energies somewhat 
below the latitude cut-off (5) is still almost correct at some not too small depth 
below the top of the atmosphere, say at Jungfraujoch or sea level. 
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We obtain for the probability of emission of m mesons with the last emission 
} act taking place in the interval dx, : 


E dE’ Ee E” dE” 
Ody, =yE? he = =|, Pia |e 0 (Fe) Ndi eee (6) 


We first carry out the ae over £, using the variable Y =log (E/E’): 


Pain BS] asa | Pu aesp{—(+))Y\a¥ | “o (F) aB"ONEe, 


‘The integration over the y,,’s can now be performed very easily using (4’): 


a Pn-1exp{—(y+1)Y}dY= Mang {{ Ym) exp {—(v + 1)yn} dy," aI : 


(n—1)! 
We denote for abbreviation 


: PRN, (Nd 
Da - \dy = a 
yt i: oly)exp {—(y + Iyjdy | o (=) (=) a 


: (eG) Got Dart (8) 


which depends on the elementary distribution w(«/E) only and is a pure number. 
‘The remaining integrations in (7) yield 


0° dE’ (Eo /E" srg Bgl dE" ao 
tat t= fem 1 
=—| of) aes i (9) 
Collecting our formulae, we obtain 
n—1 
pdx = Ga exp | ON ural, 4) DINK, so ey (10) 


This has now still to be integrated over x, because the last act of emission 
(and all previous acts) can take place anywhere between 0 and x, say, if x is the 
total length of path. The probability for no emission to take place between 
x, and x is unity because the energy is there less than F,, and therefore x, must 
not be replaced by x in the exponential exp(—@Nx,). Thus we have to integrate 
(10) over x, from 0 tox. Denoting the incomplete [-function by 


&@ gn—l As r 1 
Pu ee (11) 
we get finally CL MONe Oi ea) ce (12) 


This p is the probability for a nucleon to commit within the distance » in 
nuclear matter exactly m acts of meson emission and afterwards to have an energy 
less than £,. 

Next we consider the second case where the nucleon has, within the distance x, 
committed m acts of meson emission but has retained an energy greater than F,. 
This probability Pr, say, is obtained from (4) directly, replacing n—1 by 1, 
x, by x, and averaging over the primary spectrum (5). We obtain immediately 


n= —_ —— exp(—®Nxjot., we ee (13) 
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The total probability for emission of m mesons within the distance x is the sum: 
Pn=Pn+Pn ve eees (14) 
p® can be split up, for a more symmetrical representation, 
me = w+ i{T,(ONx) = is +(PNx)}, bia a ese ni as (15) 
and hence 
Pa=OMZal (ON) — Op Va OND). © Venn (16) 


The last step in the calculation is now to consider the variation in the length of 
path x. Ifthe nuclear diameter is d,,the probability fora path length x is 2xdx/d4. 
This is to be multiplied into (16) and integrated over x from 0 to dy. Changing 
the order of integration, the integration is elementary and gives (we denote the 
result by P,,) 


n(in+1 i 
IP = Ory] {Pulaa) oF ae Py ola} 


— Oya i ¢1(@4) — wees) Tys(a)} y 


E Vd. 
a-ONdw ora= fi (5)(1-5) F 


(17) is our final result. As can be seen, it depends on (i) the total cross section ® ; 
(ii) a weighted average over the energy spectrum of the elementary cross section, 
41, the weight function being (1—e«/E)’; (iii) the power y of the primary 
spectrum (y~ 1-5); (iv) the atomic weight, through d,. For sake of complete- 
ness we also give the formula for Py (n=0), because (17) does not immediately 
permit the transition tom=0. One gets 


2 
Po= GTiQa)- tees (17’). 


It can easily be verified that (17) and (17’) are normalized to unity: UP, =1. 
0 


§3. SOME GENERALIZATIONS 


It is highly probable that amongst the m mesons emitted some are neutrettos 
which are not observed directly in either the cloud chamber or the photographic 
plate. We shall be more interested in the probability, P} say, for the emission 
of » charged mesons than in P,,. If the probability is known for the particle 
emitted in an individual emission act to be a neutretto, 5 say, the probability for 
v charged mesons to be emitted is derived from P,, by* 


agen ey ene 
Re aa (Ty aa) aad aie pian een Gln 3t 


n= 


* We do not distinguish here between proton—proton, proton—neutron, neutron—proton etc. 


collisions. It is by no means evident, or even probable, that the cross sections for these collisions: 


are all equal and that in all of them the probabilities for the emission of a neutretto are the same 
fraction. In fact Peng and Morette (1948) arrived, for low energies, at the result that the cross 
sections are quite different, although the charge-symmetrical theory is used. Nothing is known 


about this point for the energies in which we are interested. Since, however, any such differences. 


must average out in the end, when both neutrons and protons are considered as “‘ primaries ”’,. 
because there is certainly some symmetry in the proton—-neutron system, we may assume that the: 
probabilities of neutretto emission are the same fraction of the total, irrespective of the charge of 
the nucleon. 
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The evaluation of this sum is straightforward and gives a formula for P’ that 
differs from P,, essentially in a re-definition of the constants. We get 


Prax {o'r | Toe) - 9 rai) | 


oft | Pala) eNODy Jay) got (19) 
with _ 
=ONd ,(1 —dw,,,), Digg = Vag Oty al non hn cincais ie? CLF) 
1-6 
Xy+1 = sala : 


Apart from one common factor, y,,,, which will be quite close to 1, this means 
that w is to be replaced by yw and ® by (1 —dw). 
For Pj one finds 


v7 a Mu 
Phat ayy i t2x 42 =e ay eee eee (19”) 


® is here the total cross section, including that for emission of neutrettos. In the 
charge symmetrical theory 5 has probably the value 4 and, as will be seen in § 4, 
the value of Oy 41 is between 0-5 and 1, probably § orso. x4, is then close to 1. 
Likewise a, is only slightly reduced. Since the vite of ® and w,,, are anyhow 
not accurately known, we shall work below with formula (17) and evaluate it for 
a range of values ® and w,,,. Since we shall renormalize our results to the total 
number of showers observed, the constant factor x drops out. 

A second generalization concerns the contribution of fast recoil protons. 
The experiments on the positive excess of shower particles suggest that as many 
as half the number of shower particles are protons (Rochester 1949). This is - 
not surprising. We expect indeed that one fast recoil nucleon is produced for 
each meson produced. On the average half of these will however be neutrons. 
Since neutrettos are also accompanied by recoil nucleons we expect for four charged 
mesons three fast recoil protons. This figure may be further increased since 
the strong interaction of the recoil nucleons with the nucleus gives rise not only 
to a general heating up of the nucleus but also to a few further fast ejected nucleons. 
Evidently the mechanism of the production of recoil particles is very complicated, 
much more so than for the mesons, which can be expected to leave the nucleus 
unhindered (the most that could happen to a meson is a nuclear scattering). 
Instead of using some model for the ejection of fast protons we adopt a more 
phenomenological procedure. 

Since in the experiments only the total number of fast particles is measured 
(so far), we would prefer to have the probabilities for a total of N particles, mesons 
+ protons, rather than those for a given number of mesons. Obviously a certain 
statistical correlation may exist connecting the number of mesons with the number 
of protons. (In the first place one recoil nucleon is produced for each nucleon.) 
If that correlation were quite strict, each meson would be accompanied by one 
proton, and we should merely need to divide the experimental total number of 
particles by two to obtain the number of mesons. ‘This is one way of comparing 
the experiments with the theory (cf. § 6), but it is certainly an exaggeration. 
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Alternatively we may assume that the fast protons are statistically independent 
of the mesons. We assume then that each observed particle has a probability €, 
say, of being a meson and 1 —€ of being a proton, independently of how many 
particles there are. The experiments indicate €~4. We have already the proba- 
bilities P’ for a given number of charged mesons. And weask, what are then the 
probabilities Py, for a given total of particles to be ejected? ‘This procedure may 
give too wide a spread of the number of protons relative to that of mesons. Py 
can be calculated as follows: 

Supposing PX were known, on the hypothesis of statistical independence, 
then the P;’s would be determined by 


NZ 


Las (") Ps arse y ant phew hee ae (20) 


To give the answer to our question, we have to solve (20) for the Px’s*. 
20) can be solved for the Px’s. The solution, as caa ue be verified, is 
N 2 J 


ee a CD ai Ae ora eee art) oc 21 
N am ¢ ( ) ue EH(] —£)N-# ( ) 
It is remarkable that, es for a given N the number of mesons is certainly 
less than or equal to N, P%, is expressed by the P’’s with «>N. This shows that 
there can be no arbitrariness in the P’’s. (21) has the same form as (18). We 
get then the explicit expression for Pt again by a modification of the constants: 


Prxpaaqosiat| x(a) SFr vse’) | 


ae 
WN " N+ 1 N +2 . 
—af% [Pye POO, tai) |b, eee (22) 
at — =a, te =e), 
, 1 ” Biel} , (1 —8)w 


YIN BEDE PRS ipa Caer Each BOE) Sw et 25 
dropping the index y+1 in» and y. We get for PX again the same formula, 
only with different values for the constants. Using this formula below, we shall 
Put e=—-s. 


§4. THE NUMERICAL VALUES OF © AND ay41 


Since w(e/E) is by definition normalized to 1, w,.,, must belessthan1. If we 
tentatively put w(«/£)=(B+1)(1—e/E)’ (the true law will probably be well 
approximated by this with some value of 8), convergence at « = FE requires B>—1. 
It is almost certain that B>1 because there is no known law of energy loss in 
physics that makes high energy losses more probable than small ones. Then 
+11=(B+1)/(B+y+1). As was mentioned in HJ, the theory of HHP gives 


* It is far from certain that (20) has always solutions with all P%’s positive, when the P/’s are 
given arbitrarily. In fact the assumption of statistical independence a the P»’s and the probability 
3 for being a meson or proton may be self-contradictory. In the present case, when the P;’s are 
given by a formula of the type (19) no difficulty arises. 


ofom 


¢-meg ge = 
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‘results close to B=2*. With y=1-5, we see that w,,, would be 0-4 for B=0, 
but this is also improbable. For B=2, w,,,=%. It is quite safe to say that 
-W,.+4 lies between 0-5 and 1 and, more probably, near.0-6 or 0-7. If neutrettos 
are included, and if we accept, for instance, w,,,=%, we get x,1,=+ (for 5=4 
and w/.4=%. If we calculate the total number of particles N, including protons, 
we get (with €=4) y,,,=14/11 and @) 4 = x)41)4,=8/11 =0-73. If we exclude 
memtrettos sand accept (o,.5—@,4= 4%, (we get x,.,=6/5, wi.,=7=08, If 
the low energy mesons mentioned in the footnote immediately preceding exist, 
-w, w’, ” are still slightly larger. When comparing the results with the experi- 
ments we should therefore choose some fairly small value for w, about 0-6 to 0-7, 
when only mesons are considered, and some higher value, 0-7 or 0-8, when all 
‘particles are considered. 

The total cross section © must be expected to be of the order of a few times 
‘(f/uc)?. From the absorption coefficient of fast (meson-producing) nucleons a 
‘relation between ® and w,,, was determined in HJ, namely 


O(1 —w,,,) =1-15(%/pe)? C= 2801 ae ear ace (23) 
In © neutrettos are, of course, included. Putting, for instance, w,,, = 3, we get 
® =3-5(h/y0c)?. 


This agrees, as was already mentioned in HJ, with the predictions of the HHP 
‘theory, within the limits given by that theory. 

The quantity Nd,(h/ue)? has the value 0-95 for oxygen and about 1-9 for 
Ag or Br. Thus we should get, if neutrettos are emitted, a,=3-3, a',=2-6 for 
-oxygen, and a,=6°6, a/,=5-2 for Ag or Br. For the calculation of the total 
number of particles according to (22) we have a’, =4-1 for oxygen and 8:2 for Ag, 
Br. ‘These figures increase further if the low energy mesons mentioned in the 
‘footnote below exist (as may well be the case, though not in the exaggerated 
way mentioned). We have then for oxygen a4=6-6, a, =5-2 for mesons alone, 
a’, =8:-2 for the total number of particles, and for Ag, Br twice these values. In 
all, we can say that for the calculation of mesons alone some values are to be 
-expected as follows :—a@’,~2-5—5 for O or C, ~5-10 for Ag or Br, and for the 
‘total number a’,~4~-9 for O or C, ~8-18 for Ag or Br. 

Finally we must also remark that in the derivation of the connection (23) the 
-contribution from recoil nucleons to the effective absorption coefficient has been 
neglected. It is hardly likely that this will amount to much because ® has already 
about the largest value it can possibly have, namely a few times (f/uc)?. In fact 
we have shown in HJ (in ‘‘ Note added in proof’’) that this effect is probably small. 


* Actually this theory predicts also the production of mesons with fairly low energies, of the 
order 3 x 108 ev., which do not follow the law (1). Their number is on the average equal to the 
number of fast mesons, but their contribution to the energy loss is of a smaller order of magnitude. 

“We have not considered these mesons in this paper. This could easily be done if we neglect the 
-energy loss due to these low energy mesons completely (although this is, of course, a gross exagger- 
ation) and add a term to (1) of the form ®6(¢/E) de/E. It can then easily be seen that the net effect 
is 2 slight increase of the constant w to (w+1)/2, i.e. from % to 3 (if the numbers of slow and 
-fast mesons are equal) and an increase of @ by a factor 2. The distinction between these two 
energy groups of mesons is, of course, highly artificial, and the true distribution is certainly a smooth 
transition between both groups, but the consideration shows that nothing changes, apart from a 
. change of constants, if the law (1) is not strictly obeyed, as long as the main contribution to the energy 
loss follows from (1). 
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All this serves only as a guide to the values of the constants to be expected; 
exact values cannot be predicted, but if the experiments can be explained with 
some values of w and a, similar to the above we can say that this is physically 
reasonable. 

§5. NUMERICAL RESULTS 


We give the numerical results for Py or P,, according to formula (17) for a 
variety of values a, and w,,,. The formulae (19) and (22) differ from (17) only 
by a re-definition of the constants, the constant factors y,,,, or x;,,; amount merely 
to a renormalization of all probabilities for 40 (the cases n=0 are in any case 
not observed). The range of values for the constants given, a,=2-5—20, 
w =0-6—0-9, are wide enough to include the cases considered. w,,,,is independent 
of the atomic weight. 

The probabilities are normalized to 100, including the case m=0, but Po is in 
all cases less than 20%. 

P,, decreases more or less rapidly with increasing . It is remarkable that for 
small n, P,, varies very slowly with a,, i.e. the atomic weight and the total cross 
section. It is the number of large showers which depends sensitively on ® and A. 

We see that small showers with one or two mesons are far the most frequent. 
(Here, however, our approximationsare very crude.) In light materials the shower 
size is limited to a maximum of 3—5 mesons (larger numbers being extremely 
improbable), in Ag, Br the largest showers to be expected (apart from rare 
exceptions) are 6-10. The total number of particles (including protons) is larger, 
as the larger values of constants have to be used. 


§6. COMPARISON WITH EXPERIMENTS AND DISCUSSION 


To compare our results with the data, we consider only showers with at least 
three thin tracks, for the following reason. The small showers are evidently 
produced by primaries with energies not much above those in the region where 
the cross section drops more or less rapidly to a negligible value with decreasing 
energy. Here the approximation made in § 2 is too crude. ‘The same approxi- 
mation had been made in HJ in the derivation of the power spectrum at all depths. 
The power spectrum may therefore not hold in this low energy region either 
(indeed the decreasing energy loss may make these primaries relatively more 
frequent). We cannot therefore expect that the present calculation gives the 
number of small showers very satisfactorily. 

Of the thin tracks we take it from the positive excess (Rochester 1949) that 
roughly half are protons, the rest mesons. We do not expect any of the thin 
tracks to be electrons, although they accompany as a rule the larger showers in 
cloud-chamber pictures. ‘These however are very probably the decay products 
of some short-lived type of meson which the photographic plate would not show 
to originate from the same grain (unless the lifetime were extremely short). 

We compare our results with the experimental data in two ways. We first 
assume that half the number n’ of the observed thin tracks in each shower are 
protons. We divide then the experimental numbers 7’ by 2 and, when 7 is odd, 
attribute half the number of showers to the preceding, half to the subsequent 
even number. n’/2 =n is then the number of mesons in the shower. Thus we 
assume first a very strict correlation between the number of mesons and protons, 
which is probably much too strict. The ratio } for the number of protons is, of 
course, none too accurately known, and it is not even known if this ratio is indepen- 
dent of the shower size. The figure 4 refers to comparatively small showers. 
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Table 1. Relative Frequency (°%) of Showers P,, with n Mesons* 


w=0°6 
n 
1 2 5 4 5 6 7 8 9 10 
a4 
2:5 49 Pal 5°5 1°5 0:4 0-1 0 0) 0 0 
3 48 22 8:3 2°6 0:73 0-18 0) 0 0) 0) 
4 46 DAs eee 4-3 1:5 0:46 0:13 0 0) 0 
5 45 25 13 S74, 273 0:85 0:28 0:09 0) 0 
vi 43 25 14 7:3 3:6 1:6 0:70 0:28 0:10 0 
9 42 25 14 8:0 4:3 22 1-1 0:50 0:22 0-09 
11 41 25 14 8-2 4-6 2°5 1:3 0-69 0:34 0-16 
w=0:7 
: é 
e dose 2 53 4 5 6 7 8 9 10 11-12 13-14 15-16 
3 43 23 10 4:0 1:3 0°36 0:09 O @) 0 0 @) 0 
4 40 25 13 6:3 2:6 0-93.0:30 0:09. O 0) 0 0 0 
5 SSm oe elo Onl 63-9 1-72 0206) 0-242 0:08: 0 0 0 0) 
7 34 24 16 10 Saree Soh ales OSE (Oss MO). 9X0) 0 0 
9 38) Dee ae aN 6:7 4:1 2:4 1:3. 0:67 0:33: 0:2 0 0 
ie Bye PID. Aly ala Way Bink), Pecks als7/ al) Os) Meg) 0 0 
14 Bile 2 aly “bi (ED A Ba BO Oeil SSE) OETIS) 0:25 0:07 
w=—0°8 
n 
a (lee 4 5 6 7, 8 9 10 11-12 13-14 315 
3 By) BES Bs eyohy Pile Oooh OI 0) 0) @) 0) 0 0 
4 33 25 16 86 41 1:7 0:64 0:22 0:07 0 0) 0 0 
5 30: 24-17) 11 6-0773-08 1-4. 70-560 0221 10:07°% 0 0) 0 
Uf XG. PAl “lef VP 8-4 5-3 3-1 1-7 0°84 0°39 0-3 0 0) 
9 23 Ol Ge 2 92 66 4-4 2°8 1:-7— 1-0 0:9 0:3 0 
11 ial Siam 15a, ODT On 52d 3364 2-4. 6 1:6 0:6 - 
14 DO Semel mealcl 8:9 7:0 54 4:0 3:0 2-2 2°55 1-1 0:58 
17 Sal ARS: algal 8:7 6:9 5:3 4:2 3:2 2-4 3-2 1-6 ibe 
20 Di “lg AB aal SOL o 432) Seam 2s) 3-4 1:9 1:9 
w=0°9 


i 2 3 4 5 Cy SG 9 10 i1-12°13-14 15-16 17-13 219 


3 Sil Ry ale Wea Beil tow Mee Oeilil 0) 0 0 0 0 0 0 
4 DOS 0 wii asl G08 2 Oi 2". 0°46-0°16" 0 0 0 0 0 0 
5 Dy PB Aes es) Sep 1409 2-5 ci-2, 0°50 0-20 0 0 0 0 0 
7, Velie sett tt Sal 55034 19° £:0--—0-6 0 0 0) 0 
9 144 ee 13 11. 9407-4 5-4 3:8 2:4 2°73 04 0 9 0 
iG) Domecemen 2a 109 9.087964 5-0n 13-7 4-4 16 4 0 0 
14 ee Dey NO OS 2 506 54a or 6-9) (3°98 9 0-75 0337 
iy Vie 1025949975 980 2-" 8-5 7250'6:9"6°1  5*6° 4-8. 7-8 523 393-158 1 23 
20 iLO Smee Aa ease ti Or OS ft O tontOn 7°38 O=8)) 4st 27h 30 


* or mesons plus protons, cf. text. 
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However, as our parameters w and a, are not determined accurately either, a 
constant change in this ratio will only change the parameters. For the discussion 
of large showers, however, the possibility has to be kept in mind that the fraction of 
protons may be quite different. For the second comparison with experiments we 
assume a large measure of statistical independence of mesons and protons and use 
formula (22) to calculate the relative frequencies of the total numbers of particles NV. 

The chemical composition of the plate is such that there are 1-04 light atoms 
C, O, N to one heavy atom Agor Br. We need not distinguish between C, O, N 
nor between Ag and Br, as the nuclear diameters of these atoms are not very 
different. ‘Thus we multiply the theoretical figures of Table 1 for some rather 
small value of a, by 1-04, add the figures for about twice the value of a,, and 
renormalize the total number of showers so as to agree with the actual experimental 
figures. On account of the renormalization of the theoretical P,,’s any common 
factor x,,,, as it occurs in (19) or (22) is irrelevant. 


Table 2. Experimental and Theoretical Distribution of Shower Sizes (Mesons) 


Theoretical 
Experi- = 
E mental* | a=3 and 7 a=4 and 9 a=4 and 9 
o=0:7 w=0°6 oO=—0-7 

2 89 77 83 (41) 70 (37) 

3 31:5 42 42 (19) 42 (19) 

4 i D3 22 (7:3) 24 (9:6) 

5 17 11:5 9-6 (2:5) 13°5 (3-8) 

6 9-5 5:6 4-4 (0:8) +3. (1-4) 

7 0:5 2:7 1:9 (0:2) 3-8 (0-4) 

8 0:5 1:2 0-8 (0) 1:9 (0-1) 

9 1:5 0:5 0-4 (0) 1:0 (0) 
>10 1-0 0-3 0-15 (0) 0:75 (0) 


* Results of Powell and collaborators at Bristol. 
The figures in brackets show the contribution from light atoms C, O, N. 


Since the size—frequency distribution of showers does not seem to vary much 
with altitude, we may add the results of Powell and his collaborators for Jung- 
fraujoch and balloon heights.: Similar data have also been collected in the photo- 
graphic plate by Rochester and Page. Furthermore, Fretter (1949) has obtained a 
distribution from cloud-chamber measurements. But, as Fretter himself points 
out, some particles may be lost for geometrical reasons and also, perhaps, by 
absorption in the rather thick lead plates, and, indeed, his distribution is more 
concentrated on smaller showers than the photographic plate results. We 
therefore base our comparison on the photographic plate results alone. The first 
comparison with the theory is then given in Table 2 for a few relevant values of 
the parameters. 

The agreement is seen to be very good indeed. The actual values of these 
parameters could be determined more definitely from the number of large showers 
if more statistical data were available. ‘The number of small showers (7 =2-5) 
is given well by almost any choice of the parameters. It is to be noticed in 
particular that the large showers, containing a total of 17, 18, 22 particles, are also 
well accounted for. The theory gives for the total number to be expected with 
n>7 (w=0-6, ay=4 and 9) the value 3-3, whereas a total of 3-5 has been observed. 
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The number of showers with n=1 would, if calculated in the same way, be. 
smaller than the experimental figure. As mentioned before, we believe that 
this may be due to the crude approximations made in this energy region as well 
as to a departure from the power law at lower energies. It may even be that some 
of these have nothing to do with mesons at all but are merely the effect of the 
scattering of fairly slow (i.e. ~10° ev.) protons or neutrons by protons. 

We now carry out the same comparison with experiments but calculating the 
total number of particles and assuming a large amount of statistical independence 
of meson and proton emission, thus using the constants w”, a. These can be 
compared directly with the experimental figures. We use the latter for n>3; 
this means that we are using more material on small showers than in the former 
comparison, because there only half of the showers n=3 were attributed to two 
mesons. According to § 3 we should use now larger values for w and a”, roughly 
wo” =a’ {2/(1+a)}, a” =a'(1+w’), but of course, with the assumption of statistical. 
independence, we must expect to have to use somewhat different values. 


Table 3. Experimental and Theoretical Distribution of Shower Sizes 
N =total number of particles 


Theoretical 
N Experi- a =—7 and 14 a’=9 and 20 
mental* w=0°8 a=0°8 
3 61 55 (30) 48 (27) 
4 ie 41 (21) Some) 
5 29 Sal (Gey) 30 =(15) 
6 14 295) 22 lth) 
Z 6 AST Gi) 16 (76) 
8 5 10 (3-0) 11:5 (4-6) 
9 8 6°8 (1:5) 8:2 (2-9) 
10 9 4-7 (0-7) Soh (bez) 
fl 8 
1s 1 
14 0 2:0 (0) a 3:6 (0:6) 
15-16 0 ris 
>17 3 a 1:0 (0) 3-2 (0) 


* Results of Powell and collaborators at Bristol. 


The figures in brackets show the contribution from light atoms C, O, N. 


The agreement is almost equally good, and it is certainly not justifiable to draw 
any conclusions about a correlation of proton and meson numbers. Since 
the numerical evaluation of our formula is rather laborious (partly owing to the 
lack of suitable tables of the incomplete I'-function) we have not attempted to 
secure the “‘best fit’’ by slightly varying the constants. ‘The accuracy of the 
present data makes this hardly worth while. The experimental distribution 
shows a minimum at N=7, 8 (also in Table 2 at n=4). ‘The present material 
is compatible with the assumption that this is a statistical fluctuation and this 
has to be assumed if we wish to account for the whole distribution. (If this 
minimum should turn out to be real, a special explanation must be sought, possibly 
on lines different from those suggested here.) What appears to be certain, 
however, is that the experimental distribution shows a break in the slope at N =7, 8.. 
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"This is reproduced in the theory—not quite as pronounced as in the experiments— 
and is due to the rapid falling off of the contribution from light elements so that 
only the contribution from Ag, Br is left, with its much longer tail. 

Also in this comparison there is no reason to conclude that the large showers 
of 20 particles or so are extraordinary. One shower with 31 thin tracks has been 
found by Cosyns, Dilworth, Occhialini, Schénberg* and Vermasen, and the 
probability for this would, of course, be very small for any reasonable values 
of w and a,. However, even in this case we feel it unwarranted to interpret this 
shower by a different process, for the following reasons. 

Amongst the recoil nucleons which accompany every act of meson production 
some may be sufficiently fast to produce themselves further mesons and secondary 
nucleons. These secondary processes have been neglected throughout this paper. 
The effect would resemble a mild cascade, only probably much less pronounced 
than in the well-known electron case. The recoil nucleons will also produce 
further nucleons by nucleon-nucleon scattering if their energy is not high enough 
to produce mesons. All this tends to make an already large shower still bigger, 
and the effect is naturally more pronounced in Ag and Br than in C, O, N. If 
these secondary effects give a substantial contribution, the tail end of our dis- 
tribution curves would have to be drawn out much longer, and if the effect is 
important, the values of our constants are smaller than Tables 2 and 3 suggest 
and, instead, the larger showers owe their large size to some extent to the cascade 
effect.t Pending an investigation, we feel it unjustified to interpret even the 
shower with 31 particles as a novel phenomenon. 

An indication of the existence of a cascade effect can perhaps be seen in the 
following fact. As was pointed out to us by Professor Powell, the experimeats 
show that showers with eight or more thin tracks are nearly all produced in Ag, Br 
(as follows from the number of heavy tracks) and less than one-seventh of these 
showers are produced in C, O, N (possibly even none). The figures in brackets 
in Table 3 show that, according to the present calculations, about a quarter or a 
fifth is still due to light elements. ‘This can be understood if the cascade effect is 
important; our constants would be smaller and, considering the primary effects 
only, the distribution would fall off more rapidly than in Tables 2 and 3. 
Instead, the shower size is increased considerably, and mainly so in Ag, Br by 
the cascade effects. The fact that large showers are predominantly produced 
in heavy elements only provides a strong argument for the plural character of 
the process. If the process were purely multiple, the size of the nucleus should 
be immaterial. 

The conclusions that must be drawn from these results are clear: the statistical 
data on the size-frequency distribution of penetrating showers can well be 
accounted for by the assumption of pure plural production (with the possible 
exception of the minimum). ‘This applies also to the largest showers observed. 
There is no evidence so far for the existence of any other process such as a genuine 
multiple process. ‘The values for the total cross section for meson production 
and other constants which have to be used to represent the experimental data are 
all reasonable on general physical grounds. 


* We are indebted to Dr. Schénberg for sending us this photograph. 

+ An interesting cloud-chamber picture has been described (private communication) by Lovati, 
Mura, Salvini and Tagliaferri : two particles emerging from a penetrating shower create a second 
penetrating shower in a second lead plate. ‘This means that at least one recoil nucleon has sufficient 
energy to create a second penetrating shower. It seems that this event is not infrequent. If this 
can happen in two nuclei in succession, it surely must happen within the same nucleus also. 


: 
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On the other hand, there is no proof that multiple processes do not occur in 


appreciable numbers. With a suitably invented dependence of the multiplicity 


on the primary energy it may be possible to explain the whole of the frequency 
distribution equally well. It is even more difficult to exclude the possibility 


that some of the showers, or part of the multiplicity in a given shower, are due 


to genuine multiple processes. If, for example, the average multiplicity in an 
individual nucleon-nucleon encounter were 2 or 3, the larger showers would 
be mixed, partly plural, partly multiple, and the treatment would be much the 
same as our treatment of the proton contribution. The distribution curve 
would hardly be much different from the one obtained here (with suitable 
readjustments of the constants). 

The foliowing arguments seem to point in the direction of—at least largely— 
plural production: the maximum shower size (i.e. excepting comparatively 
rare fluctuations) agrees roughly with the number of nucleons arranged along the 
diameter of an Ag or Br nucleus. If the process were purely multiple, the average 
multiplicity would have to be a similar number, which would be somewhat 
accidental. Purely multiple processes are difficult to reconcile with the large 
cross section that follows from the observed absorption coefficient, and so some 
plurality must surely exist. Finally, as was already mentioned, large showers 
with eight or more thin tracks occur nearly always in Ag, Br but rarely in C, O, N, 
and there would be no reason for this correlation were the processes all multiple; 
this correlation is even stronger than we have reason to expect. 

All this tends to show only that the processes are probably to a large extent 
plural. In order to decide whether multiple processes contribute at all, and if so 
to what extent, further experiments are needed. The best would be experiments 
in hydrogen but, failing this, observations on the average shower size in different 
elements made separately would give strong indications. In C, O, N the shower 
sizes would be roughly half the size of those in Ag, Br if the process is all plural, 
even less than half if the cascade effects contribute, but the same if the process 
is all multiple. 


ACKNOWLEDGMENTS 


We are very much indebted to Prof. Powell for communicating to us all the 
Bristol data before publication and for interesting discussions during his stay 
in Dublin; to Dr. Rochester for sending us Miss Page’s material and for advice 
regarding the probable ratio of the proton contribution; to Dr. Schonberg for 
sending us the photograph of a shower with 31 thin tracks; to Dr. Fretter for his 
manuscript on the cloud-chamber statistics, and finally to Miss M. Houston for 
computing Table 1 for us. 


REFERENCES 


Brown, R. H., CamMERINI, U., Fow er, P. H., Herrter, H., K1nc, D.'T., and PowE Lt, C. F., 
1949, Phil. Mag., 40, 862. 

Davaporta, N., and CLEMENTEL, E., 1948, Nuovo Cim., 4, 235, 298. 

FRETTER, W. B., 1949, Phys. Rev., in the press. 

HEISENBERG, W., 1949, Z, Phys., 126, 569. 

HEITLER, W., 1949, Rev. Mod. Phys., 21, 113 (compare also THOMSON, G. P., 1949, Phil. 
Mag., 40, 589). 

HEITLER, W., and JANossy, I.., 1949, Proc. Phys. Soc. A, 62, 374. 

Penc, H. W., and Morette, C.,‘1948, Proc. R. Ivish Acad. A, 51, 17, 217. 

RocHESTER, G. D., 1949, Rev. Mod. Phys., 21, 20. 

‘THomson, G. P., 1949, Phil. Mag., 40, 5809. 


684 


The Meson Intensity at the Surface of the Earth and the 
Temperature at the Production Level 


By A. DUPERIER 


Turner and Newall Research Fellow of the University of Manchester 


Communicated by P. M.S. Blackett; MS. received 30th March 1949, and in 
amended form 11th Fuly 1949 


ABSTRACT. Meson intensity has been recorded by usinga triple-coincidence arrangement 
with a lead absorber 25 cm. thick between the counters. ‘The comparison of daily mean 
values with the height of a number of pressure-levels has led us to take into consideration 
the temperature of the air layer between 200 and 100 mb. ‘The results show that this 
temperature appears to be a factor closely controlling the meson intensity at the surface 
of the earth; the intensity increases with increasing temperature at the rate of 0-12°% per °c. 

Other factors in determining the intensity are the mass of airand the height of the 100 mb. 
level over the station. The corresponding coefficients of absorption and decay prove to be 
comparable with those obtained from other measurements for mesons of the same momentum. 
Tt has also been found that the coexistence of the three effects permits explanation of the 
variability exhibited by the barometric coefficient when evaluated from short periods of 
observations. 

These results may therefore be taken as evidence that the bulk of mesons is produced 
at the height of 100 mb. 

If the positive effect of temperature is interpreted as the result of processes of decay and 
interactions with air nuclei of 7-mesons at the production level, then the coefficient 
0:12% per °c. leads to a value of 4-9 x 1078 sec. for the upper limit of the mean life of 
7-mesons. 


Si SINT RODUCT LON 
HE results of the comparison of the cosmic-ray intensity, as measured by 
means of a counter arrangement without using absorbers, with the height 
of different pressure-levels have been reported in a former paper (Duperier 
1944). For this comparison the variation of intensity from day to day was 
assumed to be given by the regression equation 


STS sR SH Te ee ee (1) 


where J is the intensity in percentages of the mean, B the barograph reading at 
the station and yp the true absorption coefficient, H the height of the pressure-level 
and yw’ the probability of meson decay per unit length of path. The partial 
correlation 7777.z, or correlation of J with H at constant B, was then obtained by 
following standard methods of multiple correlation analysis without necessity 
for reducing the cosmic-ray data beforehand. 

Since these results were published, some other observers have made similar 
studies, but the results obtained are difficult to interpret on account of the 
procedure followed. Benedetto (1946) has compared the meson intensity at the 
ground with the height of different pressure-levels, but after reducing the data 
by applying the so-called barometric coefficient 8. This is the coefficient which 
results from assuming that meson intensity is a function of atmospheric pressure 
only. It is given by B=6//5B. But if Jis in reality represented by equation (1) 


— 
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then, by substituting, we have B=y+y'5H/8B. This expression shows that 
f by being a function of 8H/8B will vary with the weather, since the factor 5H/5B 
depends on the change in atmospheric temperature accompanying the change in 
pressure 6B. It is therefore clear that the data after being reduced by applying 
f will appear to be related in one way or another to the height of a certain pressure- 
level according to the meteorological conditions during the period of observations. 

Millican and Loughridge (1948) have compared the intensity of mesons with 
the temperature of the air mass underlying different pressure-levels, but as they 
have also reduced the data first by using 6 the results are again uncertain for the 
same reason. 

It will be seen from the present paper that there is another mechanism, in 
addition to that expressed by the above expression, which also causes the 
coefficient 8 to depend on the weather. 

An account is given in this paper of a new study, similar to that made formerly, 
but using now a lead absorber of 25cm. thickness between the counter-trays. 
The comparison of the meson intensity with the height of a number of pressure- 
levels has led us to take into consideration the temperature of the air layer lying 
between 200 and 100mb. The rather unexpected result is found that this 
temperature appears to be a factor closelv controlling the meson intensity at 
ground-level; the latter intensity is found to increase with the temperature of 
the layer. 


§2. MEASUREMENTS AND METEOROLOGICAL DATA 


The meson intensity has been measured by using a counter arrangement 
consisting of three trays of counters of three units each. The effective length 
and the diameter of the tubes were about 60cm. and 5cm. respectively. The 
upper and intermediate trays were directly above each other and separated from 
the lower tray by a block of lead of 25cm. thickness. This lead was supported 
on a steel platform lcm. thick. The counting rate was about 7,400 triple 
coincidences per hour, and was recorded photographically every hour. 

The following five periods of observations have been analysed separately : 

(1) 1-4, 7-11 Jan. 1948 (9 days), 
(2) 21 Jan.—9 Feb. 1948 (20 days), 
(3) 12-29 Feb. 1948 (18 days), 

(4) 27 June-3 July, 5-24 July 1948 (27 days), 

(5) 5-10 Nov., 12 Nov.-8 Dec. 1948 (33 days). 
During the first three periods the apparatus was installed on the top floor at the 
Imperial College of Science and during the last two on the top floor at the Birkbeck 
College, London. 

For the analysis, the hourly numbers of coincidences as well as the barograph 
readings at the station have been averaged in groups of 24 hours centred for 
convenience at 0030G.m.T. The barograph was daily compared with a barometer 
of mercury. 

The upper-air data have been extracted from the Daily Weather Report of 
the Meteorological Office, London, using the radio-sonde observations made 
every day at 03, 09, 15, 21 hours at Larkhill (about 100km. to the south-west 
of London) and at Downham Market (about 75 km. to the north-east). ‘Thus the 
height of a given pressure-level which corresponds to each 24-hour group is 
the mean of eight different determinations. This procedure his been considered 
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necessary in order to minimize the meceorological errors arising from the radio- 
sonde itself and from rapidly changing weather conditions. The comparison 
of the means of the four determinations for each day at the two stations has 
shown however that these means are very closely related. Even in the case of 
the 100 mb. level the correlation, then appearing at its lowest, has a value not 
less than 0-70. 

When the radio-sonde balloons do not reach the upper levels the missing data 
have been obtained by interpolation after considering the data from the nearest 
ascent and the simultaneous ascent at the other station. ‘The error so introduced 
in the average value for the corresponding 24-hour group is negligible. No level 
higher than 100 mb. has been considered because the number of balloons reaching 
these levels is much reduced. 


§3. ANALYSIS OF THE DATA 


To begin with, the assumption has again been made that the variation of meson 
intensity from day to day is a function, and a function only, of the variation in the 
mass of air over the station and the change in height of the pressure-level at which 


Table 1 
Period 1; 7,,=—0-96 Period 2; r;,=—0-93 
Pressure- Pressure- 
level (mb.) "IH ‘BH "IHB  "IBH level (mb.) "TH "BH "IH-B ‘IB-H 
100 —0-75 0:66 —0-:52 —0-94 100 —-0°81 0-74 —0-49 —0-85 
200 —0:95 0:87 —0-84 —0-86 200 —0:97 0:89 —0-82 —0-59 © 
500 —0:97 0:92 —0-77 —0-72 500 --0:98 0:94 —0-84 —0-22 
900 —0:98 0:99 —0:65 +0-27 900 —0:94 1:00 —0:52 +0-38 
Period 3; 7,;,=—0:69 Period 4; r,;,=—0-69 
100 —0-49 0:04 —0-63 —0-76 100 —0:16 0:01 —0-:25 —0-70 
200 —0-75 0:26 —0-:81 —0-77 200 —0:49 0:17 —0:52 —0-71 
500 —0°85 0:39 —0-87 —0-73 500 —0:74 0:60 —0:56 —0-47 
900 —0-83 0-89 —0-67 +0-20 900 —0-74 0:99 —0-43 +0:27 


Period 5; r,,=—0:87 


Pressure- 


level (mb.) "IH "BH "IH-B "IBH 


100 —0°55 0:24 —0-72 —0-85 
200 —0°85 0:64 —0:77 —0-81 
500 —0:94 0-84 —0:78 —0-45 
900 —0°92 0:94 —0-60 —0-02 


mesons originate, that is to say, it is again assumed that equation (1) is valid. 
Then by comparing the partial correlation coefficients r747., and 77.7 for different 
pressure-levels it will be possible to see whether the assumption is justified and if 
so to obtain the pressure-level at which the bulk of mesons is formed. 

Table 1 gives the values which have been obtained for these correlations 
together with those for 7;,, 7; and rj, (total correlations) for each period, 


———— 
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To facilitate the comparison, the values in the last two columns of the table 
have been plotted in Figures 1 and 2. It can be seen that in spite of the low 
values of some coefficients, the change of 7;4., with the pressure level, as well as 
the change of 7;p.,4 is practically the same in all the five cases. The significance of 
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figure 1. Correlation of mesen intensity with height of different pressure-levels at constant 
ground pressure. 


In Figures 1-3 the figures near the points represent the period of observations to 
which the points correspond. 
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Figure 2. Correlation of meson intensity with ground pressure at constant height 
of pressure-level for different pressure levels. 


both changes could not therefore be questioned, but the fact that 77. increases 
continually with height while 7;,,., increases only up to about 500 mb., with a very 
pronounced decrease from 200 to 100 mb., seems to indicate that the variation 
of meson intensity at ground level is not accounted for entirely by equation (1). 
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If this had been the case we should have expected the two correlations to have — 


reached their maximum at about the same level, that for which equation (1) would 
have been valid. Another factor, therefore, must be involved. We will assume 


that this factor is the density of the air in the proximity of a given pressure-level. _ 
For a given pressure-level P we will consider the mean density of the layer 


between P and P+100mb. Instead of the actual density, it will be convenient 
to make use of the temperature of the layer, this being inversely proportional to 


the density. We therefore assume that the variation of meson intensity from day — 


to day at the ground is represented by the regression equation 
ol = p0B-wOll +a0l, 8 Sasa (2) 


where the two first terms on the right-hand side have the same meaning as before _ 


and T is the mean temperature of the layer between P, the chosen pressure-level, — 


and P+100. 

The partial correlations of J with this particular temperature which have been 
found by using the third period of observations are given as an example in Table 2. 
The Table gives also the total correlations between the different variables (with 


mee 


omission of those already given in Table 1) which are necessary to obtain ry;7.g7- _ 


Table 2 
Pressure- 
level (mb.) "BT "Hr ‘IT ‘IT-BH 
100 —0-38 —0:56 +0-80 -+0-68 
200 —0:28 —0:96 +0:74 -+0-01 


500 --0"11 --0-92 —0-67 -—0:07 


As the Table shows, the total correlations 7,7 for the two lower pressure-levels: 
are purely accidental, since they both disappear when the other two variables are 


eliminated. Instead, the correlation for the highest pressure-level remains with — 


its value only slightly reduced after the elimination of B and H. 

The same result has been obtained for the 100 mb. level with the data for the- 
other four periods. This is shown in Table 3 where the results from period 3. 
have been inserted to facilitate comparison. 


Table 3 
Retiod’ © Tay "AT "rr "IT-BH 
1 = 087 —0:56 0:92 0-75 
2 —0-:70 —0:30 0-75 Os7s 
3 —0-38 —0:56 0:80 0-68 
4 —0-25 —0:67 0:64 0:74 
5 —=0-70 —0-°58 0-88 0:63 


The fact that all the five periods give the same high value for r;7.37, makes the. 


reality of this correlation unquestionable. Even when taken separately, each 
value proves to be statistically significant. Following Fisher’s method, the 


probability that correlations of this magnitude should arise from random sampling — 


of uncorrelated observations is about 0-02 in the case of the first period, and less 
than 0-01 in the case of the others. 


In addition, this correlation allows us to explain the rather remarkable change 
with pressure-level of the correlations rz. illustrated in Figure 1. If J is closely 


related to T’ we have to expect that the correlations of J with H will be similar to 


those of T with H. In Figure 3 has been plotted the correlations 774, which have: 
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been obtained by using the upper-air data for two of the observation periods and 
for the same pressure-levels as in Figure 1. It may be noticed that the change 
with pressure-level in Figure 3 is entirely similar to the change in Figure 1. 
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Figure 3. Correlation of the temperature of the 200-100 mb. layer with the height 
of different pressure-levels. 


The conclusion must therefore be that the temperature of the air layer 
between 200 and 100 mb., or its equivalent, the density, is a controlling factor 
for the intensity of the penetrating component of cosmic rays at sea-level, this 
intensity being the greater the higher the temperature of the layer. (According 
to meteorological observations the mean height of this layer is about 14km. 
_and its average depth 4-2km. At mean and high latitudes it lies in the lower 
stratosphere; in the equatorial belt it is just below the tropopause.) 

For the partial correlation of meson intensity with the height of the 100 mb. 
level and with ground pressure, the following values have been obtained : 


Table 4 
Period "IB-HT "TH-BT 
1 —0-84 —0-74 
2 —0-58 —0-79 
3 —()-72 —0:33 
4 —0-71 +0:28 
5 —0-82 —0:49 


It can be seen that with the exception of r;4.,7 for the fourth period, all the values 
are physically of the right sign and sufficiently high to make also these two 
correlations significant. "The exceptional value +0-28 may be partly due to the 
slight variation of H during the fourth period, which could accentuate the 
importance of the meteorological errors. As shown by Table 6, it is just for this 
period that the standard deviation of H is at its lowest. In any case, however, 
if the five periods are taken together, the correlation 7;4.,7 proves to be negative 
and statistically significant, with a level of significance of 0-01. 
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When the multiple correlation 77.py7, OF seerelation between the two members. ; 


of equation (2), is considered, the following values are obtained : 


Table 5 
Period 1 » 3 4 5 


trany 0°99 ~«=—«098 = 0-9 0:88 0-96 


Equation (2) appears, therefore, to be justified when the pressure-level of 
100 mb. is taken. 


Figure 4, corresponding to the two first periods of observations, illustrates. 


the daily course of meson intensity and the three meteorological variables. 


The coefficients of the equation as well as their probable errors can be deter- © 


mined by using the relevant correlations together with the standard deviations 
o of each of the variables. These are given in Table 6, and average values have 
been obtained as follows : 
—(1-:05 + 0-16)% per cm.Hg, or —(0-77 + 0-12) x 10-3cm?/gm. 
pw’ = —(3-90 + 1-10)% per km. 
a= +(0°1234 0-024)% per °c. 


Table 6 
Period o(%) o,(cm.Hg) o,(m.)  a7(" C.) 
i Wie333) 0-62 73 4-48 
2 1-68 0-86 129 4-44 
3 0:67 0:39 63 3°82 
- 0-64 0-44 48 3-54 
5 1-75 0-66 83 4-50 


To check these coefficients in the absence of comparable determinations only 
indirect methods seem to be possible. 


Coefficient of mass absorption. 


For the coeflicient y. we can use the results of measurements of the absorption. 
curve in water. According to Ehmert (1937), the absorption coefficient in water 
down to a depth of 45 metres is given by p = 1-56/h where / is the depth in gm/cm?.. 
To use this formula in the present case the range of mesons in air and lead as. 
compared with the range in water must be considered. By using the Bethe—Bloch 
formula and applying the Fermi correction, Neher and Stever (1940) have found 
that the range in air is the same as that in water for mesons of rather more than 
10°ev. energy. From the theoretical results of Rossi and Greisen (1941) it may be 
seen that the range in lead is about 1-6 times greater than in water for mesons of 
10° ev. energy. When vertical incidence only is considered, the mass of water 
equivalent to the mean barograph reading during our observations (757-5 mm.Hg). 
together with the thickness of the lead block, steel platform, roof and counter-box, 
then appears to be 1,530 g/cm*, which would lead to a coefficient 


1-56/1,530 1-0 x 10-8 cm2/g, 


Had, however, the aperture angles of the apparatus been taken into account: 
(18° x 63° to the vertical) this value would have been a little smaller. The two 
coefficients are, therefore, quite comparable. 
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Decay coefficient. 


As mentioned before, .’ represents the probability of meson decay per unit 
length of path. The average value —3-90°/km. has been obtained after rejecting 
period 4 because the correlation of J with H was positive for this period. The 
mean range before decay would be L =(25-6 + 7:2) km., which permits evaluatiop 
of the mean life of ordinary mesons. 

‘The cut-off momentum of the apparatus corresponding to about 28 cm. of lead 
equivalent is estimated to be 45x 10%ev/c. From this value and by using 
Janossy’s expression for the average momentum of the mesons along their path 
down to sea-level (Janossy 1948), we obtain p=3-5 x 10%ev/c. As for the mass, 
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Figure 4. Variation of daily mean values of meson intensity and the three 
atmospheric variables. 


we can take some 200m., the former measurements suggesting an average of 
180 m, while more recent determinations by means both of photographic emulsions 
and cloud chamber observations give about 215m,. Thus from L=pro/m we 
have for the mean life (2-49 + 0-70) ysec., a value which is also comparable with the 
value 2-15 psec. given by Nereson and Rossi (1943) which is generally accepted. 

If the multi-directional character of the radiation recorded by the apparatus 
had been taken into account, it is clear from equation (2) that u’ would have been 
a little smaller, or L greater, as 6H increases with the secant of the zenith angle. 
The value for 7), however, would have been but little affected since the average 
momentum then to be used would have also been greater. 


Positive temperature effect. 

The possible meaning of « will be discussed in the last section. It may however 
be noted at this stage that the increase of the cosmic-ray intensity with the temper- 
ature of the 200-100 mb. layer (third term of (2)) is not incompatible with the 
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decrease of intensity with increase of the temperature of the atmosphere as a 
whole. The latter results from the rising of the meson-producing level caused by 
the thermal expansion of the atmosphere, and is represented by the second term 
of equation (2). Thus we must suppose that two separate processes are taking 
place in the 200-100 mb. layer, one of which leads to an increase and one a decrease 
of cosmic-ray intensity at sea level with rising temperature of the layer. The 
latter negative effect is almost certainly duetothe decay of u-mesons, and of course 
also occurs throughout the atmosphere down to the ground. It will be shown 
that the former positive effect is probably due to the decay of a short-lived 
intermediate product, probably in fact to 7-mesons. | 


§4. VARIABILITY OF THE BAROMETRIC COEFFICIENT 


Further support to the above results is given by the fact that equation (2) 
also provides an explanation of the variability of the so-called barometric 
coefficient. 

This is the coefficient which results from ascribing the whole variation of 
cosmic-ray intensity to the variation only of atmospheric pressure. It is given 
by the regression equation 


Si=p8B=r, SBtenon ee (3) 


It is well known that 8, when determined from observations covering only 
short periods, say 20 days, shows not infrequently a variability much greater 
than would be expected from the accuracy of the measurements. In the results 
obtained by Lindholm (1944) differences up to 100°% and more occur among the 
average values for different periods. In the same way we have great discrep- 
ancies among the values given by different observers. To quote only a few: 
3-1°%/cm.Hg with a shield of 36 cm. Pb (Barnothy and Forro), 1-1% with 20 cm.Pb 
(Messerschmidt), 1:6% with 12cm.Pb (Compton and Turner), all quoted by 
Rossi (1939); 3-4% with 20-5cm.Pb and 1-5 with 30cm.Pb (Smith and Gast, 
quoted by Mullican and Loughridge 1948). 

This variability is much reduced when long periods are considered. The 
results obtained by Hogg (1947) from a period covering more than five years 
show that the whole range of the variation of monthly mean values is only 33°% 
of the mean for the whole period. 

Our results give another example of the variability of 8. From the relevant 
data given in previous tables the following values can be directly obtained as 
indicated by expression (3) : 


Table 7 
Period B(°%/cm.Hg) ’ RH "RD 
l —(2:06-+0:22) 0-66 —0-:87 
p = (82220407) P0742" a0 
5 —(2:33+0:-37) 0:24 —0-:70 
Mean —(2:22+0:16) 
3 —(1:16-+0-30) 0:04 —0:38 


4 (1 00250:21) 0-01 =0°25 


Mean —(1:07-.0-16) 
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It may be seen that the first three values are about twice as great as the last 

‘two. Now ifthe second term of equation (2) is substituted for 5/ in (3) we have 
BSP eon eB aol, 2 Meats (4) 

which shows that 8, by depending on the variations 8H and 57 accompanying 
‘the variation of pressure 5B, will change with the meteorological conditions during 
‘the period. From the correlations given in Table 7 it can be seen that the values 
of 6H/5B and 57/5B for the first three periods would not only be relatively high 
but of the right sign to make all the three terms on the right-hand side of (4) 
negative, thus explaining the high value for 8. In fact the values of the second 
-and third terms for these periods are —0-36% and —0-61%/cm.Hg respectively. 

On the other hand, for the last two periods 5H/6B is nearly zero, and 67/5B 
‘is very small, with the result that 6 is only slightly greater than the coefficient 
of mass absorption. 

It may be noticed that had the term «57/5B been ignored it would not have 
been possible to account entirely for such a change in the barometric coefficient. 

The fact that the values of rp are all negative and those of rgy positive is in 
-accordance with the meteorological result that the stratosphere tends to be 
warmer over cyclones and colder over anticyclones, whereas most of the tropo- 
sphere is colder when the pressure is low and warmer with high pressure. On 
-account of this tendency, it is to be expected that over long periods the barometric 
coefficient will be numerically greater than the coefficient of mass absorption ; 
‘the observations confirm this. 


§5. CONCLUSIONS 

The above results show that the variation of meson intensity from day to day at 
‘the surface of the earth is due to fluctuations in air mass and in the height of the 
100 mb. level and in the temperature in its proximity. They may therefore be 
taken as further evidence that the bulk of mesons is formed at or little above this 
level. 

From a study of the results obtained in different parts of the world of the 
-effects of ground temperature on cosmic-ray intensity, we had already come to 
the same conclusion regarding the height of meson production (Duperier 1948). 
It was then suggested that the gradual increase of intensity up to the height of 27 mb. 
as found by Schein, Jesse and Wollan might be due to primary protons passing 
‘through the apparatus. Recently Pomerantz (1949), by investigating the variation 
of cosmic-ray intensity with altitude in a series of free-balloon flights by using 
different amounts of absorber between the counters, has been able to evaluate 
the relative intensity of the electron and proton-plus-meson components. He has 
found that the maximum of the proton-plus-meson curve occurs at the 80 mb. 
level while that of the electron component appears to be located a little lower at 
-about 104mb. ‘These results by Pomerantz seem to suggest that a maximum of 
‘the meson component alone must occur somewhere between these two levels, 
which should be in agreement with our result. 


So> INTERPRETATION OF THE POSITIVE EFFECT 
OF TEMPERATURE 


If, as is suggested, the penetrating component is formed by the actions of 
“primary cosmic rays on nuclei of air, the positive effect of temperature in the 
-proximity of the production level appears to indicate that the number of mesons 
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at the surface of the earth increases as the density of nuclei in that region decreases. 
This result although apparently anomalous may be explained, at least qualitatively, 
if we accept the suggestion made by Lattes, Occhialini and Powell (1947) that the 
greater number of the mesons observed at sea level are the result of decay of 
m-mesons. Following their view we shall assume that 7-mesons originate in the 
actions of primary rays on nuclei of air and that thereafter they may either be 
captured by other nuclei with which they collide or disintegrate into 4-mesons 
and presumably a neutral particle. The probability of capture would be the greater 
the shorter the mean free path of the 7-meson for a nuclear collision and the 
longer its life-time. If R is the mean free path in gm/cm? before a collision takes 
place and if L is the range before decay in centimetres, the probability of capture 
would be given by L/(R/p) where p is the density of air. If L/(R/p) is small, the 
probability of not being captured, that is the probability of decay of the 7-meson, 
would therefore be P=1—Lp/R. The temperature coefficient would therefore 
be 
éP L dp 
en ee ee ee (5) 


which agrees in sign with the experimental effect. 

For a quantitative check we can proceed to estimate the life-time of the 
m-meson and compare the results with other experiments. 

If m is the mass per unit area and h the depth of the 200-100 mb. layer, 


of 


> 


p> and sa- oo. 


From meteorological data 
f= 4-2 kms dh/oT = 20-5.mi* c., and m=102¢m. 
Consequently 
oe 2 LOC i ee (6) 


As it appears that it is the temperature of the 200-100 mb. layer only that has 
an appreciable effect on the meson intensity at sea level, it is to be assumed that 
the process of decay and capture of the heavy meson occurs mainly on this layer, 
thus indicating that the mean free path of the heavy meson is of the order of 
magnitude of the depth of this layer. We might, therefore, take R to be of 
the order of 100 gm/cm?. ; 

Owing to the very small value of the vertical gradient of temperature in the 
stratosphere it might be expected, however, that the results obtained in this paper 
would have been but little altered if a thinner layer starting down at the 100 mb. 
level had been considered. This has been tested by using the first and the second 
periods of observations for the 150-100 mb. layer; the results obtained show that 
the correlation of meson intensity with the temperature of this layer and the 
coefhicient « are both practically the same as in the case of the 200-100 mb. layer. 
In fact, they appear to be a little higher, but the difference is not statistically 
significant. ‘The value 100 gm/cm? for the mean free path of the, heavy meson 
might therefore be taken as an upper limit. 

For the average energy of the ordinary meson at the production level we may 
use 5 x 10°ev. (derived from Blackett’s spectrum at sea level and assuming a 
total energy loss in air of 2x 10%ev.). Allowing for the energy taken by the 
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neutral particle originating in the z-disintegration we have E,~=10ev. If 
finally we take for the mass 286 m, and substitute 0-12°% for « we obcain from (6) 
TX 4-9 x 10-8 sec. 

This is consistent with the results from other experiments. Camerini, 
Muirhead, Powell, and Ritson (1948) have obtained 6 x 10-sec. by comparing 
the relative number of z- and p-particles in the upward stream of mesons near 
the ground. Richardson (1948) has given, from other experiments, 1-1 x 10-8 sec. 
for the mean life of negative 7-mesons. Finally, Greisen (1948), by using the 
results of measurements at great depths underground, has claimed the value 
62 lO? sec: 

Another consideration can be made which leads to further support for this 
interpretation. According to (5) the positive effect of temperature, being pro- 
portional to L, should increase with the energy of the particles. Records of the 
cosmic radiation under 40 metres of water by Rau (1939) show that the intensity 
had a maximum in June and a minimum in January, the difference between these 
extreme values being about 4%. As the average energy under 40 metres of water 
is about five times greater than at the surface, we should expect for the temperature 
effect at this depth 5 x 0-12 =0-60%/°c. approximately. 

Now, from the Ramanathan diagram for the vertical distribution of temperature 
of the atmosphere as modified by Gregg (1931) and by Palmen (1934), the change 
in temperature at 100 mb. from winter to summer in mean latitudes can be 
estimated at some 5 or 10°c. Using Rau’s results this would give for « some 
0-8 or 0-4°% per °c., which-is quite consistent with the predicted value from the 
result at the surface. 

However, the existence of heavier mesons, as observed by different workers, 
makes it doubtful if the process of meson formation is as simple as has been 
visualized. In particular, the recent report by Brown et al. (1949) of a meson 
of 900m, which in disintegrating gives rise to 7-mesons would seem to suggest 
that this heavier particle may be partially, at any rate, responsible for the effect 
of density at the production level. 


* * * * * 


In the experiments referred toin §1, when no absorber was used so that the 
total radiation was recorded, it was found that the variations of intensity at sea 
level were well accounted for by the fluctuations of pressure and by the height of 
the 100 mb. level, that is, the effect of the temperature of the 200-100 mb. layer 
passed unnoticed. Owing to the smaller value of the average energy of the mesons 
and to the presence of the electronic component due to the lack of the absorber, 
it is possible that the greater effects of decay and absorption in this case may have 
contributed to mask the effect of temperature of the 200-100 mb. layer which, 
in addition, appears to decrease as the energy of the mesons decreases. 

Another possibility is the following. If the electrons of the soft component 
at sea level are partly the product of meson decay in the lower part of the atmo- 
sphere, then it is to be expected that the variations of intensity due to the vertical 
motion of the 100 mb. level will be partly offset by the electrons resulting from the 
disintegration of mesons. ‘The fact that these variations appear to be well 
correlated with the variations in height of the 100 mb. level seems therefore to 
indicate that something happens or tends to happen in the opposite direction 
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that compensates for the increase or decrease of the electron component. As 
mentioned above, the thermal change in the troposphere which accompanies a 
change of pressure tends to be opposite to that in the stratosphere. When the 
100 mb. level rises owing to the warming of the troposphere the number of mesons 
decreases not only because of the rise in height but also because the 200-100 mb. 
layer becomes colder. The latter decrease would compensate for the increase 
in the number of electrons, thus making the intensity appear well correlated with 
the height of the 100 mb. level. 

The former results with the total radiation might therefore be regarded as not 
inconsistent with the results as given in this papar. 
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ABSTRACT. The production of 7-mesons (mass about 900 electron masses) by cosmic - 
rays in rather thin layers of air in the upper atmosphere seems to indicate that the nucleon . 
7-meson coupling constant is not very small. The Berkeley experiments show that this is | 
also true for the nucleon—7-meson coupling constant. Consequently the decay of a charged 
7-meson into a photon and a charged z7-meson by means of virtual creation and annihilation 
of nucleons ought'to be fairly probable. In this paper this decay process, with the emission 
of one photon, is studied for the various types of r- and7z-mesons. Unless both mesons have 
zero spin, the matrix elements, and therefore the lifetimes, are found to be of the same order - 
of magnitude, and the 7-meson turns out to be too short-lived to be observed photo- - 
graphically. 

If both mesons have zero spin the emission of only one photon is forbidden. The 

detailed discussion of higher order decay processes in this case is reserved for a later paper. 


Sf INTRODUCTION 


HE decay of a charged meson into a lighter meson and a photon by means 

of virtual creation and annihilation of nucleons was discussed by Finkel-- 

stein (1947). He carried out a perturbation calculation using as the inter- 
action energy between meson and nucleon fields only that part of the total 
interaction Hamiltonian which contains the derivatives of the meson wave 
functions. Assuming both kinds of meson fields to interact reasonably strongly 
with nucleons—as is the case for the Schwinger mixture of nuclear forces—and 
taking the mass of the final meson as 177 electron masses, he found a lifetime of 
about 10-1*sec. for a vector meson at rest decaying into a pseudoscalar meson 
and a photon. 

It will be shown that if that part of the total interaction Hamiltonian not 
containing derivatives of the meson wave functions were used as interaction 
energy in the calculation referred to above, the position would be improved 
considerably. With the same parameters the lifetime then becomes 107!*sec. 
However, it is now certain that the lightest meson (u-meson) of a mass about 200 - 
electron masses is not strongly coupled with nucleons so that a 7-meson (~300 
electron masses) decaying in the way described by Finkelstein would have a very 
long lifetime, and may therefore decay in some other way. Butin connection with 
the discovery of a meson having a mass of about 900 electron masses (Powell 
et al. 1949), Oppenheimerf has pointed out that the production of these 
mesons in fairly thin layers of air leads one to expect it to have a strong coupling 
with nucleons. In that case Finkelstein’s calculations would predict a lifetime 
too short to be compatible with the observation of tracks of these particles in the 
photographic plate. 

* Now at the University of Natal. 

{+ Unpublished, remark made at Physics Conference at Birmingham in September 1948. 

Note addedin proof. ‘Through private communication the author learned from Professor Powell 
and Dr. B. Peters that at present there is no reliable evidence for the existence of short-lived 


mesons of mass around 900 electron masses. However, this does not affect the validity of this 
discussion of the possible properties of such particles. 


698 C. B. van Wyk 


The purpose of the present investigation is to see whether this result is 
dependent on the character of the mesons which are taking part in the decay process 
and on the type of the coupling. 

If both the initial and final mesons have spin zero, the emission of one photon 
.is forbidden by the conservation of angular momentum. In this case the lifetime 
will be determined by a process of higher order which will be discussed in a 
subsequent paper. 

In all other cases the difference between coupling terms with or without 
‘derivatives of the meson field variables, to which reference has been made, is less 
important here since the dependence expresses itself through a power of py/M 
where jp is the rest energy of the initial meson and M that of the nucleon. Since 
ithis ratio is about $ for 7-mesons, this factor is not very important. 

Apart from this dependence it will be shown that in all cases the lifetime is 
very short compared to 10-1” sec. so that the difficulty pointed out by Oppenheimer 
~would remain. 

Unless both 7- and z-mesons have spin zero, one would have to conclude 
either that the present formalism is not adequate, or that the generation of 
7-mesons must be due to a process quite different from that usually assumed. 


§2. NUCLEON-MESON INTERACTION TERMS 


In the perturbation calculation only interaction terms linear in the coupling 
‘constants will be used. For charged scalar mesons this term in the Hamiltonian 
-is (Kemmer 1938) 


(477)! [O* {eB + f/x[0/0xy +a V]} {tpnU + TxpU*}O dx, 


where g and f are constants with the dimension of charge and k=mc/h=p/hce; 
m=meson mass; © is the nucleon wave function; U is the meson wave function; 
-a, 8 are the Dirac matrices; tpy and Typ are operators creating a proton and a 
neutron respectively. 

Imagine the system to be enclosed in a cubic box, side length Z. Introduce 
‘the expansions 
O=[-8D eco EC ee van yo) nue ee | 
u’, (1) aN L oO) J \u"*() XP \ o*o)s |> 


l r=1, 2 


U =(heL-*)* & (2k,)-*U_,(k) exp {t(kx — Roxp)}, 
k 


where the upper line refers to the neutron component, and the lower line to the 
proton component, of the nucleon wave function and where v, u, U(k) are the 
occupation operators for neutrons, protons and mesons respectively (the signs 
refer to charge in case of protons and mesons and magnetic moment in the case 
of aneutron). a’(/) and b’(/) are spinors with the property 


{(ad) + BK jal?) = lat), 
((x)t— BK} (Z) = hb(2), 


where = +(/?+K?)?, K=M/fic, Mz=nucleon rest energy. The matrix 
‘element then becomes 


(277/L*e)*hicu*(Q)gB + tp( fie) —€ + (ap) ]}u(P), 


egy aaa 
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where p=1 for absorption of a meson, 
= —1 for emission of a meson, 
u(Q) can be either a(Q) or b(Q); «=kyhic; p=khc; ®=p?+p?. 
For the pseudoscalar case the corresponding expressions are 
(47)! [O* {gBy5 + (f/x) y50/0x9 + o— ]}{tpxU +7ypU* }Od>x 
(27/L¥e)*heU*(Q)igBys +t(f/e)[ — evs + (op) }u(P), 
where V5 = —10y He %3 O=/;t. 
The interaction term for the vector meson field is 
(4m [ne{let— Ua-+ aU) + (Vo. 1 +VEO) sy 
== || conjugate complex |txp 
where (U, U,) is the meson vector-field 
VecurlU Vo = —(0/0x,)U— V Up. 
Introduce 
Us) = (eL-9))S (Aka) {ole eUs (A) S_eUea(¥) exp filha — hy)}, 
k i=1,2 
Uo(x) =(heL*)* & (2ho)*(/ x) Us (2) exp {t(kx — Roxo)}, 
k 


mhere” —r=—kK/k, (e,,r)=0, (e,, €,)=0, (e,,e,) =(€s,€,)=L - y= —7fa, 
This leads to the matrix element 
(27/L¥e)*heu*(Q){s[(€/4)(ar) —p/e] + pfP(ar) ju(P) 
for longitudinal meson and 
(27/L¥e)!heu*(Q){g(ae) + p(f/u)Ble(xe) + 20(p,e)] }u(P) 
for transverse meson. 


The corresponding expressions for the pseudovector case are 


(4m) fossa" —75U +oU) a (f/k)(Vy =) sae d2x 
HI conjugate complex TNP 
(277/L¥e)theu*(Q)tgl(</e(o7r) — (P/#)ys] —tpR(or) a P) 
(277/L¥e)theu* (Q){g(ce) + p(f/e) Bla) —te(ce)] jae P). 
These matrix elements will be used to describe the transition 
M*(0)>m*(—p) +y(P); 
where the initial meson is at rest with energy jo, the final meson has energy ¢,, 
momentum —p/c and the photon has momentum p/c. Conservation of energy 
requires wy=e,+€2=pt+e,. It is convenient to introduce 6;=«,/u) (¢=1,2); 
thus 1 =0, + 0,. 
The matrix elements can be summarized as follows: 
Initial Meson: 
S : u*(Q){goB —foju(P), 
P, : u*(Q)tigBys—Yors}u(P), 
Vi: u*(Q)go(%eo) + foB(%eo)}u(P), 
Py: u*(Q){So(seo) — efoP(eo) ju(P). 
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Final Meson: 
S : u*(O){goB + tfo(Ho/H2)[92+ 91(ar)] ju(P), 
P, : u*(Q){igeBy5 + tfo(Ho/H2)¥ 512 + (a7) ju(P), 
vy 2 W*(Q){~go(Mo/H2)[F2(ar) + 4] +foB(ar) ju(P), 
vy 2 U*(Q){go(%es) —fo(Ho/H2)B[F2(%e2) —28,0(7,€2)] ju(P), 
Ly: u*(Q){ —g2(Ho/H2)7sL92(ar) + 1] — ifeB(or) su(P), 
Toy? U*(O){ga(7e2) +fa(t40|H4a) PLA12(7 és) + 294(0e)] uP). 
The matrix element for the emission of a photon by a fermion is (xe) where e 


is the polarization vector of the photon. In all these expressions the factor 
(27/L%e)*hce has been omitted. 


§3. CALCULATION OF THE MATRIX ELEMENT 
A typical way in which the transition by means of virtual creation and 
annihilation of nucleons can take place, is the following: 


Process Matrix Element 
M*(0)—> P*(q)+N(—@) a*(q)X0(—q) 
P+(q)> P*(q)) + (2) a*(q,) Ya(q) 
N(—4)+P*(q)>m*(—p) b*(—q)Za(q;) 


where P+, N= denotes proton or antiproton, neutron or antineutron and where 
X is one of the first group, and Z one of the second group of expressions listed 
above, and Y=(ae). Conservation of momentum requires 


q=q, +P. 
The matrix element for the three-stage process is ,/4<«7d, where 
my = Xb*(— q)Za(q1)a*(qi) Ya(q)a*(q)Xb( — q), 
Acid, = (Ho e. 2eq)(€e a a €q,)> 
e,=energy of nucleon with momentum q/c and the summation is over 


the spin states of the nucleons. 
The spinors a(q) and b(q) have the Ua ets 


B a5 (a)az¥(q) = (1/2ep){ea8 n+ (9) ae + MB 
BB —gOe*(—4) =(1)2N4) 68 w~ 2)» MB 


Introduce v=q/e,; s=M/e,, 
E=(av) + Bz, e421, 
Thus m,=4SpZ(1+£,)YA+F)X(1-£). 


Five more ways in which the decay process can take place are obtained by 
permuting the three processes involving X, Yand Z. The quantitites 7; are listed 
below: 

he cae +£,)Y(1+£)X(1—-£), 
—$SpZ(1—£,)Y(1+£)X(1-£), 

m= ~AspZ(1 +F,)Y(1-£)X(1+4), 
—$SpZ(1-—£,)Y1-—£)X(1+£), 
—4$SpZ(1+£,)Y(1-£)X(1-8), 

mg=s9pZ(1—£,)Y(1+£)X(1+2£). 


=— = EE 
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The minus signs in processes (2), (4) and (5) result from the commutation 
properties of the operators u and v. 
Notice that . 
NV, 2, Vy, 21) = —,(V, 2, —V,, — 23); 
N3(V, &, Vy, 21) =1,(—V, —2, Vy, 21); 
N4(V, &, Vy, 21) = —my(—Vv, —2, —Vy, — 23) 
N5(V, 2, Vz, 21) = —Ng( —V, — 8, —Vy, — 24). 


It is therefore only necessary to calculate (say) , and 1g. 
Define 
mass of initial meson 
~~ mass of nucleon 
and expand in powers of x 
n;=n;+x4n;+...., 
d;=1—x1d;+..... 
Then n,/d; =n, + x{1n,; + °njd,;| +.x7.... 
The contribution from these six processes is 
6 
Cys.--o= (ted) © (ld). 


The expansion of the expressions d, is carried out as follows: 


a=4—P, 
or Vv, =(2,/2)v —x24I, L=piio—vr- 
Also 2 =2{1 —2x2(Wv) +2722} 4 
_ hence €g, = M/2,=€,{1 —xa(lv) + $x%27/"[1 —7v)?].... 
Also Uh De (¢=0,1,2; O9=1; e9=p) 
hence (e,+€q, <j) 1 =(1/2e,){1 +$x2[0;+(v)] +.... 
Finally dy * =4e7(2e, — py) "(eg + €g, — €2) 2 
a Feel Ont (io)| ee. ee ete 8 (1) 
da = 4<7(2e, — Ho) (Eq eS é) 
We Me UO) ewe | kts (2) 
dg =4e7(2e, + Mo) Meg teg, +41) ¢ 
=1+4we[—1—-0,4+(W)]+..-. ne ees (3) 
din: = 4°(2e, see dCs, iat ae +e) 1 
epee ern (LU) ee win avin (4) 
d5*=4e%(e, +e), +e) We, +€q,—€2) 
=1+4hw2[—0,+0,+2(lv)]+..0. nee (5) 
dg * =4et(€g +&q, +€2) Meg +€q,—€1) 
SUetNeO I 0, 2(ID)| te nw (6) 
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Additional Processes. 


Apart from these processes there are other ways in which a charged meson 


can decay into a lighter charged meson and a photon. 


Process Matrix Element 

M*(0)> M*(—p) +P) 

M*+(—p)> P*(q)+N(— 9) a*(q)X'W(—G) fF ween (7) 

Egy NCQ) 7 aie) b*(—g)Za(m) | 

Similar to (7) but with last two processes interchanged. _—......... (8) 

0 P-(-q)+N*(q) +m*(—f) a*(QZK—q) | 

M*(0)->M*(~p) +70) | oes (9) 

M*(—p)-EES( 91) + N*(q)-0 b*(— 41) X’a(q) 

0 P-(—q)+N*(q) +m*(0) a*(q)Z'b( —q) | 

ma t( 0) MSP) CP)! elt a hea Ge eM iege eeerat (10) 

M*(0)+P(—4)+N*(q)>0 b*(—q)Xa(q) | 

0 P-(—q) +N*(q) +m*(0) a*(q)Z'b(—q) | 

M+(0)+ P-(—g)+Nt(q)> 0 b*(—q)Xa(q) Se eee (11) 

m*(0)->m*(—p) +y(p). j 

Similar to (11) but with first two processes interchanged. _—.......... (12) 
‘The expressions X’, X, Z’, Z are essentially the same as for processes (1),.... (6). 


X’ means X generalized for a moving meson and Z’ means Z specialized for a 


meson at rest. 


There are still four more processes by means of which the transition can take 


place 

OP (—g) Ng) +imt( =p) +p) ~a*(g)Z" og) 
M*(0)+P tq) +N(g)=>0 b*(—q)Xa(q) 
MO) Pia} N= 9) a*(q)Xb(—q) 
Pt+(q)+N(-9)>m*t(—p)+y(p) b*(—q)Z"a(q) 
M+(0)—> P+(q,) +N-(—4) + (0) a*(q,)X"'b( —q) 
Pa(g;) FING =2¢) =p) | b*(—q)Za(q) 
=> Pl 9) IN (Qe (2) a*(q)Zb( — qy) 
P3491) + Neg) + M*(0)—>9(p) b*(—q)X"a(q) 


where X and Z are the same as for (1)....(6), but X” and Z” are derived from the 


interaction term in the Hamiltonian by replacing 
vU by {V —(de/hc) ALU, 
vU* by {V +(ce/hc)A\U*, 


where A is the electromagnetic vector potential. 


l eat (13) 
J 


J 
eecr (15) 
| ena (16) 


‘ 


’ 
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This alteration can only affect the f-terms and the matrix elements are 


Z,, =1€(fo/H2)(xe), 
Zys =1e(fo/42)(0e), 
Zi = —&(f2/H2)(Mo/ He) B{Ax (xe) +18,0(e,7)}, 
Zi, =1e( fo/M2)Bole,es), 
Zipy = — fa/2)(Ho/ M2) B{P2%(e,7) — 18,(2e) }, 
Ziow = e(f2/ 2) Bate ,e2). 
‘To get X’’, which refers to the absorption of a meson, the only change is the 
sign and the replacement of f, by /p etc. 
The matrix elements for the individual processes are as follows: 
, =4SpX"(1-B)Z(1+E,), 
n, =zopX'(14+F)Z(1—£,), 
nj =}SpX'(1+£)Z(1-E,), 
Nyy =t5pX(1+ £)Z'(1—£), 
mj, =}SpX(1+B)Z’(1—B), 
Nig=tSpX(1—£)Z'(1+E), 
m3=_SpX(1+£)Z"(1—E), 
Ny4=5pX(1—£)Z’(1+£), 
Nys= 4SpX"(1—E(Z(1 + £)), 
nyg=}SpX"(1+E)Z(1— Fy), 


where 7/....n{2 must still be multiplied by the matrix element of the process in 


4 


which a meson emits a photon to get m,....m,). Since this matrix element 
vanishes for a scalar or pseudoscalar meson finally or initially at rest, some of 
these matrix elements need not be calculated in these special cases. 

The contribution from processes (7), (8) and (9) is | 


Cr99 = (1/2€,) {1 [(93 — 95) ][m7/d, + g/d] + (1 |2¢,)(15/dp) + ) 
where 9; =(u5+)")*/bo 93 =H2/Ho 
ae = Ze (eq =e Chae é,)* 
=1+42[0,+(Ww)]+.... 
dg * = 2e(€, eg, + Oa) 2 
=1+4re2[-—6,+(lv)]+.... 
dy t= 4eG(eg teq, +2) (eg t€q, + Fah) * 
=1+4xe[ —0,-—0,+2(lv)]+.... 


3 


Similarly 
Cyo, 11,12 =(1/2«,) {(1 |2€ q)(% 9/419) — [1/p4o(1 — 92) N33 /4ay + 2/442] } > 
where dj,!= 4e2(2e, + be) 1(Zegt+ po)? 
=1+4x2(—1-03)+.... 


=2,(2¢,+ M2) 1 =1—4w203+.... 


7 
dig = 2e(2eg— fy) Ll ther +..... 
48-2 
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16 
Similarly C13, 14,15,16 = — (1/2«,) leat 


where dig = 2e(2e, +) 1=1—gaz+.... 
dj =2,(2e,—po) 1 =l+dez4+..-- 
dyg) =2,(e, +<g, —€2) 1 = 1+ 4x2[0,+()] +..-- 
dig! = 2e,(e, +e, to) t= 14 $x2[ — 0, +(lv)] +... - 
The reciprocial of the lifetime of the meson is given by 


q1=(2n/h)| A|% ‘ 
where 


H= (27/ L?)3/2(he L)8 | [(<p€1€2)*(27hc)?] faratc, oeeeg + C3 ee a £2) + Cis ocee opt: 


and p=density of final states per unit energy. 


§4. ADEQUACY OF THE LEADING TERM 


The series obtained for 4(n,/d;) by expanding the numerators and denominators. — 


in a power series of x is certainly convergent. 
It will be shown that only the leading term is important, even if x=4. 
Notice that 
vi(V, z)= = Vil = Vie), 


2,(v, 2) = —2,(—v, —2). 
Anticipate the form 
n=(z/2)A +(1 + 2,/2)B+x2,C, 


which will be derived in the next section. Suppose the expressions X and Z 
from which 7 is derived contain only one type of Dirac matrix like S(«p) and not 
a linear combination like B(ap)+AP(ck). Then A and B contain the same 
vectors v,e;.... and therefore (z,/z)A and (1 +2,/z)B behave similarly under a 
change in sign of vand z. From the above this implies also a change in sign of 
v, and z,. Because of the additional factor (7) in the spur represented by C, this 
expression contains one v more (or less) than either A or B. 

Hence n(v, z)= +n(—v, —2). From the series for 2, follows that in the 
expansion 

n=n+x1n+.... 


consecutive terms behave differently towards a change in sign of v. The same is 
true for dy*+d,'. If n,(v,z)=n,(—v, —3) then 


n,(v,2)= —m,(—v, —2)=—n,(v, 3), 


hence n,d>'+,d,'=n,(d>!—dy') which, after integration over the directions 
of v, results in an odd or even series of x. Closer investigation shows that the 


second term in the series for nd, '+m,d,' is smaller than the leading term by about — 


a factor $x22707< 4-32? (G4 Oe <1): 

The same argument applies to the groups (2, 3), (5,6), (7, 8,9), (10, 11, 12), 
(13, 14,15,16). This result holds also for those linear combinations of different 
types of matrices used in this paper. 


On the Decay of t-Mesons 705 


Example. 


Before calculating the matrix element of a particular decay process it is 
convenient to eliminate v, from the calculation. The equation 


V, =(2,/2)v —x2,1 


is invariant under the various changes in sign of v,z,v, and 2, so that after 


elimination of v, the calculation of m, and ng is still sufficient to determine 1, 13, 
N,N; also. 


(421) + Bay =(24/3){(av) + Bz} — x2, (a/). 
Hence 1+£, =(z,/z)1+£)+1—-(z,/z) —x2;,(a!). 
Since Y=(ae) and (1+ £)(«e)(1 + £) =2(ev)(1 + £), therefore 
my = 3(ev)(2/2)SpZ(1 +E)X(1—B) 
+ 3(1—2,/3)SpZ(ae)(1 + E)X(1—E) 
— $x2z,SpZ(al)(ae)(1+ E)X(1— £). 
Since (1 —E)(ae)(1 + £) =2(ae)(1 + £)—2(ev)(1 +), 
therefore Me, = — t(ev)(2,/2)SpZ(1+F)X(14+£) 
+4(1+2,/2)SpZ(we)(1 +E)X(1 +2) 
+ hxz,SpZ(al)(ae)(1 + £)X(1+4+ £). 
It is evident, from the relationship between n; (¢=1....4), that 3 N= 254 (V, 2) 
where /(—v, —2)=f(v, 2)... Similarly n; +n, = F(z, v, 2), where am 
(2, —V, ~2)= — Fe, Vv, 2)4 
Consider now the transition 
My (0)>m,(—p) +y(P). 
A vector meson at rest decays into a pseudoscalar meson, momentum —p/c and 


a photon with momentum p/c. 
To calculate the term in egpg. put 


X = (xe), Z=1By;. 
SpZ(1+£)X(1 —£) =0, 
SpZ(ae)(1 + E)X(1 — £)=0, 
SpZ(al)(ae)(1 + £)X(1 — £) = — 2izSpy,(al)(ae)(aeg) 


= 82(e,ép)l. 
4 
‘Thus 2, = — x22,(e,é9)/ leading to — & n;=0. 
i=l 
4 
However X n,fd,=n,{d;'+dz1-—d;'—d,*}=3x2n,. 
i=1 


Similarly n,=O=n;. ‘Thus neglecting higher powers of x 
6 
x n,/d;= — 3x7(e,e9)r0,2°. 
i=l 


The ef, fs-term was calculated by Finkelstein and is of order x° (see Table). 
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It can easily be seen that the additional processes do not contribute. 
mn, =1SpX'(1-£)Z(1+£,), 


where X’ =(ae}), Z=iBy;, if in the intermediate state the heavy meson is a trans- 
verse one. Because of the y; in Z 


n+ = const. (€9,v)v =0 =ng = 5, 


Similarly if the heavy meson is longitudinal in the intermediate state. 

To get m9, 1,4, the expressions {9,711,212 must be multiplied by a factor 
which vanishes if the secondary meson is pseudoscalar. 

Processes (13)....(16) only contribute to f-terms. Thus 


AH =(27/L? poe 1¢ 9)? X?(€y,€)70,(M. [Ar )ego£, 
p= 450, L?/(2rhe)®, 


giving 
Th = {(eqne)r}?{x°0M | (87°6,h) }{e°gog3/(tic)? 


Summing over thetwo directions for e and averaging over the three directions of €p, 
71 = x30) Merg?93/{127?h8,(hic)?} =1012 


for the parameters used by Finkelstein. 

The picture presented here is however less attractive if it is required to describe 
the decay of mesons of mass about 900 electron masses. 

To get an idea of the order of magnitude in such a case, put x =4, 6, =4,0,= 2, 
giving 7 1=6 x 10 gip5/(fc)?. The lifetime r= 10-°sec. would therefore require 
gig3|(hc)?=10-8. Since the Berkeley experiments seem to confirm that g3/Zic is 
not very small and probably about 1/10, it would follow that g3/hc would have 
to be less than 10-7, which seems incompatible with a reasonable cross section 
for creation of r-mesons unless some quite new process comes into play. 


Table. 


The matrix elements for all choices of initial and final mesons are tabulated 
below. The scalar, pseudoscalar, vector and pseudovector mesons are denoted 
by s, ps, v,pv, respectively. Longitudinal and transverse vector mesons are 
distinguished by lv and ty respectively, similarly for pseudovector mesons. 
The symbol (7, 7) denotes the group of four matrix elements containing g,g;, 
8b jf: Sif; (im this order) in the decay process where 7 is the initial, and j the final 
meson. ‘T'he second column contains the total matrix element for the first six 
processes, while the contribution from the rest appears in the third column. 
The entry q’ in this column denotes a quadratic divergence. The expressions 
containing a factor z* are finite when integrated over q, but those containing only 
zg are logarithmic divergent. 

For the sake of consistency all divergences will be ignored. This prescription 
will in some cases require the calculation of new leading terms. 
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d d 
Process RO te [Cree 
(s, tpv) | —(z/3)x?(e,e.)1{3 +46, —100,2?}23 i 
2x(enee)l(Hl}4a){1 —(1-+0,)22}2 
2x(e,es)l2° 2x(Ep,€o)10o3° 
41(e,e2)l(Wo/[2)% 
(s, tv) 4(ee,)2° 
2x(2€2)(Ho/t2) {81+ (1 —40, + 63) 2? } 2 
(21/3) x(ee2)(A2—8,)(z" + 2)z q 
4i(e€s)(€o/49)z 41(€e2)(Uo/[42)(82—9,)2 
(ps, tpv) | —4z(ee.)z 
— 2x(ee2)(Mo/2){8, + (1 — 6, +65) 2? jz 
(27/3) x(ee,){2(1 —26,) + (1+70,)2?}2 q" 
4(e€2)(€2/p2)2° —4 (ees) (Uo/H2) 072° 
(ps, tv) 3x7(E,en)12° 
2x(Epe2)1(Uo/[42){1 C= 0,)27}z 
—A4x(e,esylz* 2x(€,€2)1052° 
—4 (en €2)l(Mo/2) 2" 4 (ep €o)l(€o/p42)2° 
(v, s) 4(ee9)z° 2(e€))ab?6,2° 
— (21/3)2(eeo)(Ho/[42){2 + (1 + 6?) 2?+ 26424} 2 q 
2x(e€)){ —8,+(1+20,)2?}z x(ee,)ab?0,2° 
— 41 (ee) (€2/[42)2 — 41 (ee )(Mo/P2)2 
(v, ps) —3x?(e,eq)lz* 
2x(€p€0)l(Mo/[42)(2 + 82)2* 
— 2x(e,eo)l(1+ 27) 2x(ep ep )lO,2° 
4 (eno) l(po/H42)2* —4(eneo)l(€n/2)2° * 
(v,tv) | (4/3)x(eql)(ee2){2 +22? —24}z gq 
4 (eol)(€€2)(Mo/-2)% 2(€ol )(€€2){ — (Ho/p-2)a(Ge =6,)+b}z 
3x°(€ol)(€€a)2° 2(€ol) (ees) {a(O{ — 0203) + b—2}2 
(2/3)xx(eol)(€ea)(Ho/H2) {1 + (6—6,)2? + 20,24}2 gq? 
(v, lv) | (2/3)x(eeo)61(Ho/H42){4 + (3 —8,)2? + 20,24} z g° 
3x?(€€)0,2° 2(eeo) {a(97 — 6263) +b—2(u5/M3)} 61 
4(€€)0;(Uo/[12)2 2(eeo) {a(8, — 92)(4o/t42) +b}0,2 
(2/3)x(ee9)0,{1 +72? — 224 be ¢ 
v, tpv) | 22x(eo,e)eo(1+6,)2* 1(€,€2)l(€ol)ab 2° 
4 (one )€2(€2/M2)2 A (eon e)eo([o/[2)% 
4i(eonees2® 2i(€n2o)l (ol )ab3z° 
(2/3)x(Eql)(e 2nel(Ho/l42) ie ee eee 
(2/3)x(€o,e)22(€2/H2){2—8, + (1+46;)27}2 g 
v, Ipv) | 2ix(enes)r(44o/}+2)(1-+62)2" 
A(e,eo)rz 4(e,€o)r(Ho€ 2/43) 
41(E,eo)1(€o/L9) 2° —4i(eneo)rO5(Mo/42)2* 
(2/3)x(e,eo)r{2—0,+ (14 46)2?}2 he 


* This contribution cannot be neglected as Finkelstein suggested. Its inclusion lengthens the 
‘time by a factor (€4/f19)-?= 10. 
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(43) 9 

Process 43 | Gicens sep | Cr. 

(pv, s) (1/3)x?(e,eo)l{ —(3 + 46,) + 100, 27}2° i 
os 2x(epeolOs(Ho/H2)2° : 
2x(E,eo)l(—1+ 22") 2 2x(€p,€o)10 22° e 
4i(e,e0)1(Mo/ 2) 

(pr, ps) | 4t(eeo)z 2i(ee)(O;/A5)z 
— (2i/3)2(2eo)(po/t42){2 + (1 — 90, +62)x? + 26)24} 2 : ry 
2x(€€ 9) {8, — O02? kz —2x( ce) abi 2° 1 
4(€eo)(€o/[2)2° 2(e€o)(€2/p1-2)aBtz* 

(pv, tv) 2ix(€,€)€o(20, — 1)? —ix(ep€s)l(eql ab az* 

a 4i(epeo)e2(€2/42)2° 2i(epls)l(€ol)(€o/[42)2* 
4( ep €o)e22 {—2/65(e,eo)l(Eol) + 4(eneoe242}2 
(2/3)ac(eol)(epea)l(Mo/Ha) {1 + 42° —224}2+ q? 

(2/3) x(e,eo)€2(€2/42){2—8, + (1 —50;)2? a 

(pv, lv) | 2tx(e,eo)r(p4o/H2)(1 —30, + 87) 2° 
A1(e,eo)7° 
—4(eneo)r(€o/[42)2 —4(e,eo)r(Ho/2)(82—91)2 : 
(2/3)2(e,eo)r{—2+0,+(—1+56;)27}2 g? . 

(pv, tpv) | (2/3)x(e22)(eol) {4 — 522 —224}z gq . 
—4i(ee2)(€ol)(Mo/42) 3" — 2i(ees) (ol) (Uo/pW2)ab7z° 
— (i/3)x?(eeg)(eol) {3 +40, — 108,27} 2° 2i(€25)(€ol) {ab 3x” — (2/3)(1 — 22) }z 
(2/3) 2e(eee)(€ol)(H49/42) {1 — (3 + 01) 27+ 20224} g° 

(pv, lpv) | (2/3)x(eeo)(€1/m2){3 — 8, 527+ 2624} g° 
+ (i/3)x?(eep)0, {3 +46, —100, 2? }23 — 2i(eey)ab?z® 
41(ee)0,2° 2i(ee)ab}(U9/W2)2* 

+ (2/3)x(ee)0,{1 — 2222242 @ 


The notation is as follows: 
pi=meson rest energy, where the suffixes 0 and 2 refer to the initial and final mesons 
respectively. 
€.—energy of final meson. 
p=e,=energy of photon. 
e=meson polarization vector (i=0, 2). 
e=photon polarization vector. 
r=unit propagation vector of photon. 
91=P/Mo, 92=€ 2/40, 92=2/[H0, 63=(uo+p?)4/po. 
l=6>r: 
M=rest energy of nucleon. 
€g=(M?-+-q*)t. 
q/c=momentum of nucleon in intermediate state. 
X= po/M. 
2z=M/eq. 
a={83(0,—8)}*, b= {0,11—63)}4; d°q=q? dq dw. 


$5. CONCLUSION 
From the table it is clear that for x =4 all the matrix elements are roughly of 
the same order of magnitude. Recalling the result of the example it must be 


concluded that the decay process r++ 7+ + makes the t-meson too short-lived 
to agree even approximately with observation. 


ee 
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The situation is however quite different if both the initial and final mesons 
have zero spin. In this case the emission of a positron-electron pair or of two 
photons must be considered. If, in particular, one of the two mesons is scalar 
and the other pseudoscalar, then the production of a pair is forbidden. A calcu- 
lation, which will be published in a later paper, of the probability of emitting two 
photons leads in this case to a lifetime of the order of § x 10-1? sec. where F is the 
reciprocal of the product of the dimensionless coupling parameters of the two 
meson fields. This is no longer in direct disagreement with experimental 
-evidence. 

Hence the conclusion is that either r-mesons and z-mesons have spin zero, 
or the nucleon—r-meson interaction is different from that represented by the 
scheme 

P+>N 47+, 


-or the usual perturbation theory is not adequate to describe these processes. 
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ABSTRACT. The electronic states of naphthalene are determined using antisymmetrized 
molecular orbitals, which are based on wave functions given by the simple molecular orbital 
theory. ‘This method allows the electronic interaction term to be included explicitly in the 
Hamiltonian for the system. A clear distinction is made between electronic configurations 
determined with these wave functions and the actual states of the molecule: states arise 
from a superposition of configurations of the same symmetry. ‘The ground state becomes 
considerably stabilized by this configurational interaction, to which even highly excited 
levels contribute. Calculated energy values for some of the lower transitions agree well 
with ultra-violet spectral data, as do the newly derived oscillator strengths. A tentative 
assignment is made regarding the polarization of the electronic transitions in the region of 
30,000—55,000 cm™. 


§1. INTRODUCTION 

HERE are two main theories known for the theoretical discussion of molecular 
structure—the valence bond treatment (Heitler, London, Pauling) and the 
method of molecular orbitals (Hund, Mulliken). Reasonably good results 
can be obtained with either method for ground state phenomena; but for the 
interpretation of ultra-violet spectra—involving electronically excited states in 
relation to the ground state—both treatments are inadequate. Particularly 
in the polyacenes the two methods lead to different conclusions, neither of which 
can be correlated satisfactorily with experimental data. Because of this disagree- 
ment it seemed worth while to re-examine the method of molecular orbitals and 

to extend it; this has been done with reference to the naphthalene molecule. 


§2. METHOD OF MOLECULAR ORBITALS 

Our discussion of the molecular orbital (m.o.) theory falls into three stages: 
(a) the simplest form in which all electronic repulsion is averaged yields the basic 
m.o. wave functions and levels; (4) inclusion of the electron repulsion term in the 
potential energy expression leads to individual configurations; (c) configurational 
interaction gives the states of the molecule. We reserve the term state for the 
final description of the molecule, whereas any allotment of electrons to the available 
m.o. levels used in (0) is referred to as a configuration. 

The first stage deals with the method as developed chiefly by Lennard-Jones 
and Coulson (1939), who applied it to hydrocarbons. The fundamental assumption 
is that electronic interaction can be averaged over the whole molecule to give 
a self-consistent field. This means that the effective Hamiltonian contains 
one-electron functions only and may be written W = ,H(v), where H(v) is the 
effective Hamiltonian for electron v. Such an approximation enables us to use 
one-electron wave functions ¢, and to associate with each of them a molecular 
orbital energy «;. The energy of a configuration of the whole molecule is then 
simply the sum of the energies ¢; of the occupied m.o. levels. For convenience 
only z-electrons are usually considered. 


* For parts I and II see C. A. Coulson, Proc. Phys. Soc., 1948, 60, 257, and H. C. Longuet- 
Higgins, Proc. Phys. Soc., 1948, 60, 270. 
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Coulson (1948) has published detailed calculations for the polyacenes. As 
these will be required in our later work (stages (b) and (c) above) it will be necessary 
to begin by summarizing his results. In these molecules the nuclei are arranged 
symmetrically with respect to the short and long axes. Calling these the y and z 
axes respectively, we can group the nuclei into sets each of which is symmetrical 
or antisymmetrical to reflection in the planes perpendicular to the plane of the 
molecule through the y and z axes. 


BA 


Figure 1. 


M.o. wave functions are built up as linear combinations of atomic orbitals. 
at each carbon nucleus ¢;=»,c’s,, where ¢;=wave function of the 7th m.o., 
ys, =atomic 2pz wave function at carbon atom 7, (c!)? = weight with which #, enters 
the zth m.o. These wave functions themselves must bear the symmetry 
characteristics of the sets of carbon atoms.. Thus four types of wave functions 
can be written down, which we denote by S,S,, S,A,, A,S,, 4,A,, where for 
example S, indicates symmetrical behaviour with respect to reflection in the plane 
through the z axis, perpendicular to the plane of the molecule, and A, indicates 
antisymmetrical behaviour with respect to this reflection. 

In terms of the effective Hamiltonian H(v) and the normalized m.o. wave 
functions ¢; the energy of a m.o. level is determined by 


c= [O*()H(o)4i(v)dz,] [4*()b(0)4r,. 

Minimizing the energy with respect to the coefficients c’ leads to the secular 
determinant. Because of the symmetry grouping, the secular determinant 
factorizes into blocks of these symmetry types called P,Q, R,S by Coulson. 
The m.o. energies which are derived as the roots of this determinant are therefore 
denoted by P,P.P;, 0,02, R,R,R3, S,S,. Coulson gives a diagram of energy 
levels which were calculated including the overlap integral of the 2pz atomic 
orbitals at adjacent carbon atoms. ‘This diagram shows also relations between 
the energy groups P, QO, R, S, and both the symmetry character S,S,, etc. and their 
standard group theory description. 

E, is a Coulomb integral considered constant for all carbon atoms. 

In the ground state it is assumed that the five lowest orbitals are each completely 
filled with 7-electrons (each orbital then accommodates two z-electrons according 
to the Pauli principle); it therefore gives rise to a description P? Of R7 Pz Sj. 

Similarly the “‘ first excited configuration ’’—or, more accurately, the configura- 
tion which from this diagram appears to be the first excited one and which arises 
from the elevation of one electron in m.o. S, to the lowest unoccupied one, O.— 
can be described by P? QO? R? P3 S, Oo. 


Ee F. Facobs 


To determine the symmetry type to which such a configuration belongs one 
multiplies the symmetries of the electrons forming it; for the ground state 
having all m.o.’s doubly filled this gives a totally symmetric configuration. In the 
‘first excited configuration” the symmetry is given by S, x Q, i.e. A,A, x S.Ay, 
i.e. A, A complete table of symmetry types of products of two-electron 
symmetries is given in Coulson’s paper in the more usual notation of group theory 
first introduced into molecular descriptions by Mulliken (1933). 


P O R S 


ie Ajg By Boy Byg 
Q By Aig Bog Boy 
R Buy 3g Aig By 
S Byg Bou Buy Aig 


An electronic transition from the ground state is spectrally allowed if the 
product of the electron symmetries in the excited state behaves like a y or 2 vector 


, 3 2) 
' i] 


Lo 


Bog Big 


Figure 2. 


(symmetry A, or A, respectively), the direction in which the electric moment 
vector changes during the excitation, which must be in the plane of the molecule. 
Above, we deduced ¥ for the first transition; Coulson’s results for higher trans- 
itions are given in Table 1 together with values for excitation energies. 


Table 1. Naphthalene Energies (Coulson’s Results) 


Symmetry of excited configuration. By, By Bsg By Aig Bou Aig 
Excitation energy in terms of y 1:267 1532 1°715° 1:868° 1:884 2-133.5°2-30e 
Polarization y zg F z F y F 


F=forbidden. 


In this ‘Table y is a resonance integral between neighbouring carbon atoms 
and has a value of 3-4 ev. 
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In all polyacenes the lowest transition appears to be of A,,>B,, type, 
i.e. polarized in the y direction. This result is in direct contrast to valence bond 
calculations (Craig, in course of publication) which predict z polarization for the 
first allowed band in each case. 

One limitation of this m.o. method is that no account is taken of electron 
spin so that the calculated energies of excited configurations are really averages 
of the corresponding energies for singlet and triplet configurations. The 
averaging of electron repulsion uniformly over the whole molecule is a further 
approximation that cannot be justified. Particularly for some excited levels 
the calculations give electronic charge distributions that are not at all even, so 
that the potential field for each electron is certainly not equivalent at all atoms of 
the carbon framework. 


§3. EXPLICIT INTRODUCTION. OF ELECTRON REPULSION 


Our first step towards a refinement of this simple method is to include 
electronic repulsion explicitly in the Hamiltonian of the system with the m.o. 
wave functions obtained above as basis. This necessitates introducing anti-- 
symmetrized determinantal wave functions, as was done by Goeppert-Mayer 
and Sklar (1938) for benzene. 

Let the normalized wave functions for the simple m.o. be 6;=%,.c'ys,.. Then 
the wave function for the configuration is, in the first approximation, ‘’; = NI,¢;,. 
where z is taken over occupied m.o.’s and N is a normalizing constant. 

Instead of these, the wave functions of the type used by Goeppert-Mayer 
and Sklar for a configuration 7 of the molecule are 


$i(1) a(1) 
$1(2) B(2) 


b(n) Bn) 


where 7 electrons occupy m.o.’s 1,...., 7. «(1), B(1) define the spin function 
of electron 1 as «, 8 respectively, and N is a normalizing constant. This can be 
written in the form 

Vj =N2Zp(— 1)°T1,P$(u)a(e)6(r)B(), 
where P is the permutation operator permuting electrons over all available- 
da, $8 and wu, v represent electrons. 

Including the term %, ,e?/r,,(e =electronic charge, 7,, = distance between 
electrons , v)in the Hamiltonian in their investigation of benzene, these authors 
obtain fair agreement with experiment. However, all integrals except those for 
nearest neighbours are neglected and three-or-more centre integrals are omitted 
altogether (some numerical errors also occur in their integrals (Griffing 1947, 
Parr and Crawford 1948)). London (1945) has carried out similar calculations for 
diphenyl and concludes that some of the three-centre integrals may be quite 
important; they would certainly destroy the good agreement with spectral data for 
benzene (but see Parr and Crawford). Goeppert-Mayer and Sklar (1938) treat 
configurations calculated with these antisymmetrized wave functions ‘’; as giving 
the energies of the actual states of the molecule and compare the resulting energy 
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differences with those observed in the spectrum. ‘The only interactions between 
configurations which they consider arise from the special degeneracies of benzene 
that have no parallel in the higher -acenes. However, we find this to be an 


inadequate description of the energies, and further interactions between some of Y 


the calculated configurations can arise in the following way. 


$4. CONFIGURATIONAL INTERACTION 


Let the wave function ’; define the electronic configuration j of the molecule ; 
then the symmetry type of the configuration can be determined from the product 
of the symmetries of the m.o. wave functions ¢ which form ¥’;. Since several 
such configurations may belong to the same symmetry class (e.g. configurations 
_in which the m.o.’s are each completely filled are all totally symmetric) the 
possibility of interaction between them must be conceded. Configurations 
belonging to different symmetry classes or of different multiplicity cannot interact 
with one another. Parr and Crawford (1948) find in their calculations for 

-ethylene that the configuration in which both z-electrons are excited gives rise 
-to considerable interaction with the lowest configuration, stabilizing it by about 
T-SeVe 

This concept of configurational interaction finds confirmation in some recent 
work carried out in this laboratory on the hydrogen molecule (Coulson and 
Fischer 1949). For this molecule more accurate values can be obtained as fewer 
assumptions and approximations have to be made in energy computations. As 
for the z-electrons in ethylene, thereare only two totally symmetric configurations. 
It is found that even at equilibrium internuclear distance some degree of mixing 
takes place, which increases for greater nuclear separations. Only after taking 
configurational interaction into account does the ground state appear as the lowest 

. State of all, with the singly excited state falling below the doubly excited one for 
all internuclear separations. 

The fact that configurational interaction exists indicates that the anti- 
symmetrized wave functions ‘’; are not by themselves suitable for describing 
the molecule. If we consider a linear combination of ‘’;’s including a 
sufficiently large number of terms we should approximate to the true wave function 
of a state (Longuet-Higgins 1948). Our procedure is therefore to carry out 
energy calculations using antisymmetrized wave functions and to obtain the 
excited configurations and their symmetries. We then allow configurational 
interaction to take place among the lower levels in the same symmetry class. 
This third step in the calculations we consider essential as it leads to a sequence 
of energy levels differing significantly from that obtained in stage two of our method. 
Only the lower levels in each symmetry class are considered for convenience in 
computation, although theoretically all Jevels in one class should be included. 


§5. METHOD AND RESULTS FOR NAPHTHALENE 
The true Hamiltonian for the system of z-electrons is 
10 10 e 


HEHE ST) SS 


Dab v=1 v<j Ty 


w 


where 7(v) is the kinetic energy term for electron v and A(v) is the potential 
energy of electron v ina field of molecular framework when each C atom is stripped 
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of its z-electron. This A(v) differs from that introduced in § 1 (first stage of 
m.o. theory): the earlier H(v) included a term arising from electronic averaging, 
not found in the présent H(v). 

10 


We take H(v) == H,,(v) =a sum of contributions from each carbon atom k 
k=1 ; 


which has lost its z-electron. This may be written 
. (A) 2 
Ho) = Hi) —e | HET ae, 
vA 


where H7(v) is the field acting on electron v due to neutral carbon atom k. We 
further assume that the potential field around each carbon atom is the same 
(so that the form of Hf is independent of k)—an approximation here, but 
true in the case of benzene, where Goeppert-Mayer and Sklar introduced this 
substitution. 

The wave function for the ground configuration is 


Vo = (10!) #2 p(— 1)? Pgi(1)ac(1)$1(2)B(2) . . .. $5(9)x(9)65(10)8(10), 
where P, ¢, «, 6 are defined as above, and its energy is 
| Vx HWY dr, 


with corresponding expressions for other energies and with functions which are 
—allreal. It is found that all energy integrals can be given in terms of three series 
_ of expressions, «, y, 5, defined by 


e;= energy of zth molecular orbital level 


= | SFONTO) + HOM bdr, 
rae = | f ARIE) brdbolWwarsdr, 


b= |] BOW) E dolodblwendr, 


Let us use "E?? to denote an energy level of multiplicity + for which m electrons 
have been raised and the m.o. levels concerned with the excitation are pq. 
n=( will represent the ground configuration, n=1 singly excited configurations 
which alone are usually considered in m.o. calculations. But doubly excited 
configurations (7 = 2) are also important as they have the symmetry of the ground 
state and so may give rise to configurational interaction with it. Using the 
simple diagram of Coulson’s paper (Figure 2) and labelling the m.o. levels in 
order of magnitude, we expect the following excitations to be significant for 
low energy differences. 

The relation between the notation in Figure 2 and the present enumeration 
Pecueroy veveleis 7-90 5. Ky, 2 5,2.54;,0>, R, are equal to 1; 2,.3; 4, 5, 6, 7 
respectively; pP>q=5->6, 4>6, 3-6, 5-7, 4+7,....all of which have been 
evaluated below. 
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In terms of the molecular integrals «, y, 5 the energies can be expressed as 
Ty = 2(ey +g +g teg +5) + (11 + Y22 +33 + Yaa + V5) 
+4(yr0 +713 +Vi1a +15 +23 +Y24 + 05 + ¥3a +35 + 45) 
—2(8y9 +843 +844 +545 +593 +894 +595 +534 +535 +945); 
1B 96 = 2(ey +eg teg tq) +(€5 +e 6) + (11 +22 + ¥33 + ¥44 + 56) 
+4(yy2 +713 + Via +03 + ¥24 +34) +2(Y15 +16 +25 + Y26 
+ Y35 + Y36 + Yas + Y46) + 556 7 2(842 tis 513 ay O14 a 893 zi So4 “ts 534) 
— (845 +846 +595 + 5a¢ +535 +536 +545 + 5a6)> 
1 ESS = 2(ey +g teg teg tee) + (¥11 +¥02 +733 + Yaa + Yes) 
+4(yy2+V13+Y1a t+ V16 + Yo3 + Yea +Y26 +'¥34 + Ya6 +46) 
— 2849 +843 +544 +546 +523 +894 +596 +534 +536 + O4¢)- 
Similar expressions for other energies are easily written down, and in particular 
for the singlet triplet separations one obtains, for example, 
1f 56 —3F 56 — 25... 


Changing now to integrals over atomic orbitals, we neglect three-or-more 
centre integrals and introduce the notation of Goeppert-Mayer and Sklar. 


A=(1r,17)= [ $,(H,2)—ty( Wear, 
B=(0r,$8) = | bes) dW), 
C= (05485) = | $2) =H Qar, 

D= (rs 13) = | deal) =D 2)er 
Q=(H;, 5) = — | ds(AHG(Iy(1)ar, 


R=(Hy,7s) = — | ghe(1)Hs(1)y(Dar. 


Let W,,, =energy of a 7-electron in a neutral carbon atom so that 


[P(v) + Hi) }i(r) = Wophj(); 
and let the molecular orbital wave functions ¢ be expanded as 
$= & Cp}, by =X,d py, 
in which 7%, is the hydrogen-like 2pz orbital at carbon atom 7, and (c,)?, (d,)? are 
the weights with which it enters the m.o.’s ¢,, ¢, respectively. With this notation 
the following expressions are obtained : 
a Axed; a 2B. ea, Ue cya) at 22C (C74, of c.d,)(¢,d,, ati cds) ats 42D, C,CAj Ay 
oy = AXc*d? aD 2B .C,CA,d, af ZC aCe oF Cd, )(C,dy ste c.d,) te UD, (Cd, te Cyd as 
= Wop a 2x(R,.s 1p Crsler€s + x(Q ane B, (ce aR Gs): 


xTS 
At this stage the following approximations and assumptions are made: 


ip titties ee hi 


el i it lil 
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In the atomic 2pz functions which are of the form 7 sin 6e~7”? (in the usual 
notation) the screening constant Z is taken as 3:18 throughout. The C—C bond- 
length is taken as 1:39 a. everywhere, an approximation allowing the naphthalene 
molecule to consist of two hexagons with the dimensions of benzene. In the 
evaluation of the integrals C, D, QO, R only contributions from neighbouring 
atoms 7, s are taken into account, those for non-nearest neighbours being 
neglected, whereas integrals of type B (i.e. simple Coulomb terms for the atomic 
charge clouds on two distant atoms) are taken for all pairs of atoms as B decreases 
but slowly with distance between 7 and s. Corrected values for these integrals 
are used as given by Parr and Crawford (1948), with some extrapolation in type B 
integrals. Possible variations of the screening constant Z and the C—C distance 


Table 2 

(1) (2) (3) (4) (5) (1) (2) (3) (4) (5) 

0-0 - 44,, Grouhd’state| 1B2° 67876 

0:9840 ®Byy T LE,58 -6-8617 

11015 *Byy T 68683 "Aig F 
SE5® 3-1687 1E5® 6-93.49 

23-2219 8Bay, T 1B46  7-0614 
8F,57 3-4648 1E,57 7-1400 

35210 ®Byy T 1B47 7-252 

3-5624 Ay T 7-2643 Bon 4 

3-6214 By Z LE 18 7-448 

3-6503 1Boy Y 1B 177-5898 
8E2® 33-6660 7:9080 ®Boy T 
SEA? 3-7441 85309 *Byy T 
1B, 4-0657 8-6198 1Ayg F 
sE46 4.5572 2B, 8-8152 88152 *Bag T 
8E18  4-9513 1E.37 —9-1620 

49873 Aig F 1B,48 9-641 
SEI? 5-2010 | 9-9329 1Ayg F 
8B37  5-2936 10-6308 1 By, Z 

5-3972 *Ary T 10°8683) 2Bgees ) 47 
157 6 -4304 1E26 12-3728 12-3728 1Byg F 

66883 'Byy or 13-2497 1A, F 


(1) Description of configuration; (2) Energy (ev.) of configuration: zero at —4-0657; 
(3) Energy (ev.) of state : interaction taken into account; (4) Symmetry; (5) Spectral 
activity. 7, triplet; 2, z-polarized; Y, y-polarized; F, forbidden. 


in excited states are neglected. In the wave functions overlap up to two C—C 
bond distances is included. 

Using the wave functions ¢ obtained by Coulson’s method as basis for the 
antisymmetrized wave functions ’, we determine configurational energies which 
are recorded in column (2) of Table 2. Finally, resonance among configurations 
of the same symmetry is allowed for as described below, the results being shown 
in column (3). 

There are six configurations (which are shown in column (1)) of 14,, type 
whose wave functions can be denoted by ’,.,,.,. To find the energies of 
the corresponding states of the molecule we take linear combinations of these, 


PROC. PHYS. SOC. LXII, II—A 49 


718 F. facobs 


and minimize the energy calculated from the x with respect to the coefficients c,. 
Such configurational interaction should in principle take place in any symmetry 
class for all configurations within it. We have here considered it for the six 
lowest configurations of 'A4,, symmetry and for the three lowest ones in each of 
the following classes respectively: 1By,, ?Byy, 1Boy *Bou ?Ay, (two levels only). 
(‘Lower’ or ‘higher’ means here as given by calculations of method (a), § 2.) 
In the 14,, configurations this resonance is most pronounced, and it threefore 
seemed best to include levels up to about 6 ev. above the ground one, assuming 
that the influence of any higher energies on the ground state would be negligible. 
In the remaining symmetry types we assume that interactions of the lowest three 
configurations in each case will suffice as an indication of the effect. But care 
must be taken to include a sufficient number of configurations in each symmetry 
group: the differences between taking two or three levels into consideration 
(in 1B,, and !B,,, classes) are considerable both as regards energy separations and 
calculated oscillator strengths. The values in column (3) must therefore not 
be taken as final, for inclusion of still higher levels will lower these energies a little 
further. In so far as these results are valid, therefore, the calculated energy levels 
must correspond to observed bands at somewhat lower energies. 


§6. DISCUSSION 


The energy values in column (2) (i.e. of the separate configurations) differ 
from those found by Coulson both in magnitude and in the order in which they 
arise because of the corrected Hamiltonian that has included electronic repulsion 
explicitly. But this improved Hamiltonian leads to the unlikely result that it 
is energetically sometimes easier to raise two electrons than one, for the 1£,°° 
level just precedes the level 125%. That levels like 1#,°*, 1#,°* should fall fairly 
close is inherent in our procedure. We start with the simple (Coulson’s) 
molecular orbitals which give levels £,5°, F,°* well separated; we then combine 
the wave functions ¢ with spin factors «, B which separate the singly excited 
level into a triplet and a singlet level as in the diagram, bringing 1£,°* near to1F,*%, 
which itself is unaffected by spin. 


BO ee CON ee Se TCG 
Ea EE, 
1 enh 8G 
E. 56 ee 
a 7756 
Ey 


Only at the third stage—resonance among configurations of the same symmetry 
—do we reach the final arrangement of states. This resonance is strongest among 
1A,, levels and is particularly effective in lowering the ground level (by about 
4-1 ev.). Since spectral term values depend on energy differences with the ground 
state it is essential to perform calculations like those that lead to column (3) before 
a comparison with experiment can be attempted. 

In Table 3 we give a list of weights (c,”) to indicate the degree of mixing of 
configurations in the lowest two states of '4,,, 1B,y, ‘By, symmetry respectively. 
In this table configurational wave functions are listed in the order in which the 
corresponding levels occur on an energy scale (column (2), Table 2). The most 
important conclusion we reach is that this order is not always an indication of 
the importance of a level in its resonance effect. Judging by the 14,, results it 
is approximately true that the lowest configuration (1£)) is most important in the 


— 


ee 
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ground state, while the second and third configurations are the chief contributors 
to the first excited 14,, state. On the other hand, in1B,, and 1B,, states the order 
of configurations is not even approximately in accord with the weights of the 
wave functions belonging to the final states. 


Table 3 
Weights of Configurations in States 
7Aj,g configurational Ee yee Wee ao ee reas 


wave function 


c;2 in ground state 0-668 0-087 0-017. 0-063 0-155 0-009 
cp” in ist excited state 0-001 0-474 0-480 0-009 0-032 0-004 


EB configurational | y. 57 yas 16 1B, configurational hy 56 yt ar ye 7 
wave function f a wave function . , : 

‘Cy? in. 1st state 0-133 0:499 0-368 Cy” in 1st state 0-096 0-485 0-419 

.¢,” in 2nd state 0°829 0-018 0-153 c,2 in 2nd state 0-899 0-031 0-070 


Returning to Table 2 column (2), we observe that the ground level configuration 
‘does not appear as that of lowest energy, for several triplet configurations fall 
below it: apparently the triplets form a cluster of configurations preceding a 
group of singlet ones. ‘This may arise in part from too wide a singlet-triplet 
separation met with in m.o. calculations; the placing of triplets relative to 
singlets is, however, considered less significant than the order of singlet states 
alone. It is, indeed, this sequence of states which is perhaps the most unexpected 
feature of our results. In Coulson’s method the sequence-of transitions to 
singlet levels is in the following order’ 


TE PY) TEM Z) CER) TH(Z) CEPR) 1EM(Y), 


whereas now the first of all singlet levels is z-polarized, very closely followed 
by a y-polarized one. The next transition is again long-axis polarized. For 
‘comparison we quote the results of valence bond calculations obtained by Craig 
where the lowest transitions in terms of polarizations are (i) forbidden, (ii) z-polar- 
ized, (ill) y-polarized. Although this agrees with our determination that a 
g-polarized transition occurs energetically below a y-polarized one, the two pre- 
dictions do not necessarily support one another. According to our present 
calculations the first two bands are so close that they could not be detected easily 
as two distinct bands in the spectrum. Our second 1B), state (i.e. the third 
allowed transition, again z-polarized) should be attributed to the second observed 
band. This new sequence of singlet levels then gives fairly satisfactory agreement 
with Kasha’s experimental results (in course of publication). He interprets 
the spectrum of naphthalene as consisting of the 33,100 cm! weakly allowed 
transition preceded by a separate forbidden one at 31,060 cm™!. These are 
followed (in order of increasing excitation energy) by the strongly allowed 
45,000 cm! band; a further weakly allowed band appears at 52,500 cm-. 
Numerically our results for the lowest allowed transition agree fairly closely with 
the observed values, experiment giving about 4 ev. while we calculate 3-6 ev. 
Energy values for higher excited states are less reliable but should be attributed 
to bands at energies rather lower than the calculated ones. The uncertainty 
about the lowest forbidden transition persists, however; several triplet states 
49-2 
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fall immediately below our lowest allowed band, whereas the first excited 14,, state: 
(to which such a forbidden transition might well have been attributed) occurs 
1-4 ev. above the first allowed transition. 

A further test of the value of our present extension of the simple molecular 
orbital treatment is afforded by a comparison of observed and computed oscillator 
strengths. The intensity of absorption depends on the transition moment Q 
associated with that transition, and Mulliken’s (Mulliken and Rieke 1941) 
theoretical formula for determining oscillator strengths is used: 


Fneor, = 1°085 x 104vQ?, 


where v is the frequency of absorption in cm“, Q is in cm. The transition 
moment Q is calculated from the integral 


O= | x#(2r)xxdr- 


Here r is the displacement vector of the electron; y,, y;, are the wave functions 
of the states /, k between which the excitation takes place. 
A comparison with experimental values is given below. 


Table 4. Oscillator Strengths 
0-52 1:109 1-306 
(vy) (2) &) 
Vops.X 1078 (em72) 33. 45 «52 veneer. X 1072 (cem™) 29 30-5459 
02226036" 1-49 70-035 
bbeoele @ 0) @ ©) 

In this Table (a), (c) refer to the first and last of the three stages of m.o. theory 
described in § 2, so that f(a) is calculated from the simple m.o. wave functions ¢ 
using frequencies of the observed spectrum, while f(c) is determined from wave 
functions x belonging to the states of the molecule with frequencies corresponding 
to calculated energy separations. 

These new oscillator strengths represent some improvemnt over those of the 
simple form of m.o. theory. Although in the m.o. method the calculations for 
higher states are less reliable than for lower excited states, the new f values do 
reproduce the great difference in intensities of the 33,000, 45,000, 52,000 cm-? 
bands in the naphthalene spectrum. 

A further interpretation emerges from the calculated oscillator strengths f(c) ; 
the 33,000cm™ band might be attributed to a z- and a y-polarized transition nearly 
overlapping one another. This conclusion differs from results of both the 
simple m.o. theory and of valence bond theory, but it is not altogether in disagree- 
ment with experiment. In 1944 Prikhotko measured the absorption spectrum of 
naphthalene crystals at low temperatures in polarized light and found both 
z- and y-polarized lines present in the region 29,000-33,000 cm. Nor does a 
vibrational analysis of this spectral region exclude the possibility of two distinct 
sequences of lines which might belong to two separate electronic transitions. 
This dual character in the vibrational lines persists even on substitution by methyl 
groups. 


Sobs. 0-18 1-7 0:20? | ftheor, (a) 


To summarize: we consider the value of the present investigations to be two- 
fold. Firstly, they show how unreliable the results of the simple m.o. method 
are when applied to problems involving excited states; at the same time they 
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point the way to an improvement which lies in taking electronic interaction into 
account explicitly together with configurational interaction. However, the 
calculations become rather long and the new method is probably too laborious 
to be of use for more complicated molecules. Secondly, in so far as these new 
calculations have a physical meaning, they suggest a new interpretation of the 
ultra-violet spectrum of naphthalene. The ‘first’ allowed observed band may 
consist of two bands of different polarity almost coinciding with one another. 
‘The ‘second’ strong band (having its peak around 45,000 cm~) is then to be 
correlated with the next z-polarized transition. ‘The polarity of this band thus 
agrees with the results of the simple m.o. method, but it is now attributed to the 
next higher state of that symmetry. For the lowest forbidden transition the 
symmetry remains uncertain from these calculations. 
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ABSTRACT. Up to the present spectroscopic studies of flames supported on a tube 
have been limited by the geometry of the flame. Thus the important reaction zone presents: 
a three-dimensional problem in which it is impossible to locate the emitters exactly. 
Moreover this zone is so thin in premixed gas flames that intermediate products cannot be 
determined from absorption studies. "To overcome these difficulties the authors have 
constructed a burner with a two-dimensional diffusion flame in which the reaction zone at 
ordinary pressures is 5-10 mm. thick and has an optical depth of 5 cm. or more. With 
this flame many of the reacting molecules can be located through either their emission or 
absorption spectra or both. The value of the new technique in the elucidation of com- 
bustion processes is demonstrated with reference to NH,—O, and H,—O, flames, and mention 
is made of its application to hydrocarbon flames. 


$i. INTRODUCTION 


HE application of spectroscopy to combustion has been widely practised 

and has provided much of the existing evidence on the subject. It is well 

known, however, that there are a number of serious limitations in the 
methods employed, particularly from the structure of the flame. In a Bunsen 
flame for example, the most interesting reactions occur in the narrow inner cone 
and not in the large volume of interconal gas above it. The thickness of the inner 
cone, however, is usually less than 0-1 mm. and it is extremely difficult to study 
the detail. Moreover, it presents a three-dimensional problem in which it is 
virtually impossible to decide on the precise location of the emitters. The optical 
depth of the inner cone is insufficient for absorption studies and this too is a great 
drawback because absorption spectra can give quantitative as well as qualitative 
information on the molecules and atoms involved. Spectroscopic evidence 
is thus limited to overall observations such as the emission of C, and CH bands 
in the inner cone of hydrocarbon flames or to the emission and absorption of OH 
radicals in the interconal gases. 

Some of these difficulties have been partially overcome by Gaydon and 
Wolfhard (1948) who used premixed flat flames. Normally a premixed gas flame 
will strike back when the mass flow is reduced below a critical value and this occurs 
before the cone has flattened out. Under carefully controlled conditions, with 
small gas flows which approached the limiting value for self propagation, Wolfhard 
(1943) found it possible to obtain a flat flame above the burner at any pressure. 
At 1 atmosphere, however, these flames are very small and require a burner of 
about 1 mm. diameter. At low pressures (2-20 mm. Hg) the burner diameter 
may be increased to about 50 mm. and the reaction zone has a thickness of more 
than 10 mm. depending on the mixture strength. Emission spectra can be 
observed in different parts of the flame and thus the relative positions of the various 
emitters may be obtained. At reduced pressures the flat flames are too thin 
optically for absorption work and it has only been possible to detect OH radicals. 
in absorption using very high resolving power. 
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Egerton and Powling (1948) have recently developed a large flat premixed 
gas flame at atmospheric pressure by having a uniform gas flow across the burner, 
thereby eliminating the usual parabolic velocity distribution which leads to 
conical flames. The inner cone of the flame is normally very thin but widens for 
very slow and rich flames. ‘The method is very promising especially for the 
detection of radicals in the zones before and after the main reaction. 

Apart from the last example it may be said that neither premixed gas nor 
diffusion flames examined hitherto allow full use of absorption spectroscopy. The 
primary requirements for this are a broad reaction zone with good optical depth 
for absorption work. ‘These requirements have been met at normal pressure 
in a diffusion flame of unusual design which is described in this paper. 

It should perhaps be pointed out that the chemical reactions in diffusion 
flames may be different from those in a premixed gas flame. Nevertheless an 
attempt to understand these reactions is fully justified by the great technical. 
importance of diffusion flames. Apart from calculations on the height and shape 
of such flames by Burke and Schumann (1928) little fundamental work has been 
done hitherto. 


§2- DESCRIPTION OF NEW TECHNIQUE 
Diffusion flames burning with oxygen have a reaction zone up to 1 cm, in 
thickness depending on the height of the flame. This zone can be obtained in 
a more convenient form by using a flat burner in lieu of the usual cylindrical pipe. 
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Figure 2. End view of flame. 
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Figure 1. Burner. 


Figure 1 shows the design of the burner. The fuel gas and oxygen issue with 
equal velocities of about 20 cm/sec. from parallel rectangular tubes (5 cm. x 0-7. cm.) 
which have one long side in common. The gases rise vertically and diffuse 
into each other in a laminar flow on account of the low Reynolds number. 
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On ignition burning takes place at the interface and a flat diffusion flame is 
obtained which presents a two-dimensional problem, rather similar to the premixed 
flat flames described in §1. The flame front is kept very steady by stabilizing 
it with a surrounding stream of nitrogen which flows up the outer jacket. ‘The 
optical depth of the flame is the length of the interface and is arbitrarily fixed by 
the burner dimensions; in the investigation described here it was 5 cm. This 
depth at ordinary pressure is sufficient for absorption studies. Quartz windows 
are suitably sited in the outer jacket. 

This burner has not been used so far for a detailed study of combustion 
problems, but an assessment of its potential usefulness has been made by finding 
out what bands can be detected in absorption and emission and the location of 
these bands in the reaction zone. 

Combustion occurs where the fuel and oxygen are in contact at the interface. 
Figure 2 shows an end view of a H,—O, and an NH;-O, flame. The reaction zone 
is thin at first but broadens higher up the flame. ‘The combustion products, 
viz. H,O or H,O and N,, etc., are formed in the reaction zone and diffuse outwards 
on both sides whereas the oxygen and the hydrogen (or ammonia) diffuse towards 
the reaction zone and maintain combustion. The reaction zone of the hydrogen 
flame emits in the visible region a fairly uniform blue continuum similar to that 
experienced with a premixed H,—O, flame, whereas the reaction zone of the 
NH,-O, flame is clearly divided into two parts. On the ammonia side a narrow 
plane of yellow radiation can beseen, which is due to the ammonia«-bands ; moving 
across the flame towards the right this is followed by a dark space and then by a 
fairly thick plane of blue continuum on the oxygen side. The flames are not 
exactly vertical even when the velocities of both gases are equal, because the 
diffusion rates of the fuel and oxygen are different. The position of the gas 
interface in the burner has therefore no relation to any point in the reaction zone 
above the burner. 

A spectroscopic study in emission was made by projecting an image of the 
end view of the flame (12 mm. above the burner and corresponding to the line 
AA’ in Figure 2), on to the slit of a spectrograph by means of a quartz-fluorite 
achromatic lens, the aperture being chosen so that the different layers in the flame 
did not overlap appreciably. For absorption work an image of a carbon arc or a 
hydrogen discharge lamp was focused on to the flame. ‘This image and the flame 
were then focused on the slit of the spectrograph. Since the most information 
about the reactants would be derived in the horizontal direction AA’ across the 
interface and not in a vertical direction parallel to the interface, the images on the 
spectrograph were made to fall across the slit either by the use of a 90° quartz 
prism or by putting the spectrograph on its side. ‘Two hairs were placed across 
the image of the flame on the slit, and cast line shadows over the whole wavelength 
region of the spectra. ‘These lines served to locate the source of the bands in 
the flame. A Hilger medium quartz spectrograph was used for most of the 
investigation. 


§3. THE HYDROGEN-OXYGEN FLAME 


The use of a H,—O, flame in this assessment is limited because hydrogen 
cannot be detected in absorption with a quartz spectrograph. The reaction is less 
complicated than the hydrocarbon reaction, or even the ammonia flame, and it 
seemed a simple application for the new technique. The results are shown 
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graphically in Figure 3. The hydrogen diffuses from the left and the oxygen from 
the right. The abscissa represents the horizontal distance (in millimetres) in 
‘the flame 12 mm. above the burner mouth. The ordinate on the right is a rough 
qualitative measure of the strength of the bands estimated from the negatives. 
No attempt has yet been made to measure the intensities accurately for calculating 
the absolute concentration of the molecules. The absorption bands of the OH 
‘radical extend over approximately 9 mm. and the emission bands are symmetrically 
located over the greater part of this distance. ‘The ordinate on the left is a 
‘temperature scale. ‘Temperature measurements in the reaction zone were made 
by the Na line-reversal method. A thin wire coated with sodium chloride was 
inserted into the flame below the point of measurement. This coloured only 
‘the small and relevant part of the flame and eliminated risk of self-reversal. ‘The 
anode of a carbon arc was used as light source together with a rotating sector 
‘to reduce the effective temperature. The results are shown by the top line of 
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Figure 3. H,—O, flat diffusion flame. 


Figure 3. The temperature attained a rather flat maximum at 2,800°c. This 
-agrees very closely with the calculated value for a stoichiometric mixture of 
H, and O,, allowing for dissociation but not for external heat losses. It appears 
therefore that diffusion flames are quite as hot as premixed flames. 

Comparing the temperatureand OH measurements, the latter cease in emission 
_at about 2,100°c. (though this depends to some extent on the sensitivity of the 
measurements). The OH absorption extends slightly further on both sides. 
By using the OH bands for temperature measurements by reversal, instead of 
‘the Na lines, the values obtained agreed within the limits of error with the sodium 
reversal temperature, and proved that in the H,—O, diffusion flame the electronic 
-excitation of molecules and atoms is thermal even in the reaction zone. This 
is believed to be the case in the premixed flame of H,-O, also but certainly not 
in premixed flames of hydrocarbons with oxygen. 

Oxygen is normally transparent in the ultra-violet, but hot oxygen absorbs 
‘because higher vibrational levels are excited which permit a transition to the 
‘upper electronic state above 2,000a. At 1,000°c. the absorption has already 
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spread to 2400a. and is readily observed. The spectrum contains a vast number 
of lines which do not form definite heads, and a satisfactory analysis has not 
been achieved up to the present. The existence of these lines however provides 
a means of detecting the oxygen in the flame. On the outside of the flame 
where the temperature is only moderate the absorption by oxygen is a complicated 
function of the temperature, but in the reaction zone where the temperature is 
above 2,000° c. the absorption must give at least a qualitative indication as to where 
the O, is used up. These results for the H,—O, flame are also shown in Figure 3.. 
Just beyond the point of maximum temperature the O, concentration had 
fallen to a value where absorption could no longer be detected. 

The height of these excited vibrational levels is interesting. The O, absorp- 
tion bands in the flame can be followed up to about 26004. ‘This corresponds 
to a vibrational energy of about 33 kcal. The reaction H +O,-OH +O which 
can be regarded as an important step in the H,—-O, reaction is normally 
endothermic (16 kcal.), but if this high vibrational energy of the O, is available it 
may be exothermic. At least one metastable O, level (+A,) is lower than 33 kcal. 
and may therefore be partially excited and play an important role in combustion. 
It seems to be a characteristic of the diffusion flame that the reactants only meet 
after they have nearly reached the theoretical flame temperature in contrast to 
the premixed flames where they are continuously in contact and begin to react 
at the ignition point (800-1,000°c.). The oxygen molecules emit as well as 
absorb, and the emission is shifted towards higher temperature compared with 
the absorption. Figure 4 (see Plate) shows an actual spectrum which gives in 
emission the (0,1) OH band, and the (0,14) and (0,15) O, bands. It is 
noticeable that both O, bands and OH bands do not occur at the same position 
in the flame. The existence of the Schumann—Runge O, bands in emission in the 
flame was rather unexpected since it is very difficult to obtain them in a discharge 
tube because the potential energy curves of the excited and unexcited states are 
relatively displaced. In a diffusion flame where the O, has high vibrational 
levels in the ground state it is easier to excite the upper state in accordance with 
the Frank—Condon principle. 

The HO, radical could not be detected either in absorption or in emission. 
The radical is supposed to play an important part in combustion and Minkoff 
(1947) has calculated that-its spectrum should lie around 2900. It has not 
been reported so far. 

In some reaction schemes H,O, has been suggested as an intermediate product 
and Egerton and Minkoff (1947) found appreciable amounts in an H,—O, flame 
at low pressure. H,O, has a continuous absorption which is not very strong in 
the quartz ultra-violet and does not make a very sensitive test. Since no 
continuum was detected down to 2100 4., however, it may be said that the H 205 
concentration was less than 0:2°% of the total pressure. 

By a similar reasoniny the partial pressure of ozone must be smaller than 
0-006°% since ozone absorbs aut ae aa near 2500 a., but this continuum was 
not observed. 


§4. THE AMMONIA-OXYGEN FLAME 
The experimental arrangement for the ammonia flame was in every way 
similar to that of the H,—O, flame. The ammonia system has an advantage over 
hydrogen in that the fuel can be observed in absorption. At room temperature: 
NH, begins to absorb at 2254. and after a few diffuse bands the absorption 
becomes a continuum towards the vacuum ultra-violet. The effect of temperature 
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on this absorption was studied separately in a quartz cell 5 cm. long which was 
filled with ammonia gas and inserted in a furnace. The length of the optical 
path was thus similar to that of the flame. As the temperature was raised new 
band heads appeared and at 1,000°c. these bands had spread up to 2440 ., the 
whole absorption becoming more continuous in appearance. Figure 5 (see 
Plate) shows how the absorption by ammonia in the cell varied with temperatures 
between 20° c. and 800° c., and Figure 6 shows the ammonia absorption observed 
in the flat diffusion flame. As the ammonia approaches the reaction zone the 
bands near 2200 a. become stronger and later continuous, then they suddenly 
cease. Figure 7 shows again in graphical form the results for the NH,—O, flame 
of all the absorption and emission spectra observed in a horizontal traverse of 
the flame 12 mm. above the burner mouth. The height of the curves is again 
only a qualitative indication of the strength of the bands. The ammonia diffuses - 
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Figure 7. NH;—O, flat diffusion flame. 


from left to right and the oxygen in the opposite direction. ‘The ammonia is 
apparently unaffected until about the 1-0 mm. point on the abscissa then the bands 
increase in strength due to the increasing temperatures and reach a maximum 
at 3-2mm. At 48mm. the NH, disappears completely. At the point of 
maximum decline in the NH, absorption, however, the ammonia «-bands appear in 
emission and these bands are usually attributed to NH, radicals. ‘The maximum 
of these bands is at 3-°3mm. At 4:0 mm. the maximum of the NH bands 
appears in absorption and the emission bands of this are strongest at 4-2 mm., 
the shift being due to the higher temperature in this direction. In this way 
the whole course of the dehydrogenation of the NH, in the flame can be followed. 

NH has been detected in absorption by Frank and Reichard (1936) who 
decomposed pure NH, at 2,000°c. NH, could not be found in absorption possibly 
because of the small dispersion of the quartz spectrograph in the visible region, 
and possibly also the multi-line nature of the band renders it difficult to find it in 
absorption. These bands are discussed again later. ‘The NH band at 3360 a. 
is much easier to find because the Q branch forms a strong head and the band 
looks like a strong line. 
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O, was again detected in absorption and spreads from the right (Figure 7) 
up to 5-4mm. Since the ammonia spreading from the left has completely 
disappeared at 4:8 mm. (see above) it appears that the ammonia and oxygen never 
come in contact and the whole reaction goes via the intermediate products. The 
relevant positions of the two bands can be seen in Figure 6. 

O, in emission is again displaced relative to the absorption on account of the 
temperature gradient in the flame. 

OH is the most prominent feature in the spectra, either emission or absorption. 
It commences near 3:6 mm. and goes far over to the oxygen side (Figure 7). 

The temperature of the ammonia flame is plotted in Figure 7 (with reference 
to the ordinate on the left) from sodium line-reversal measurements. The 
maximum temperature was found to be very close to the theoretical temperature 
for a stoichiometric mixture. The fall on either side of the maximum is similar. 
‘OH reversal temperature measurements agreed with the Na line results and it 
was again concluded that the electronic excitation of the OH and Na is purely 
thermal in the NH,-O, flame as in the H,—O, diffusion flame. Even the NH 
bands may be used to obtain the reversal temperature in the position where they 
appear and the measurement made with them was 2,200° c. at 4-1 mm. (Figure 7) 
in agreement with the other values. Figure 8 shows the position of the NH band 
compared with OH on the plate. 

Returning to the problem of the NH, bands, the following rather sensitive 
test was applied. The NH, radiation is the only strongly visible radiation and 
occurs in a narrow vertical plane less than 1 mm. thick; by holding a piece of 
white paper across the end of the flame and looking at the outside of the paper 
one should see an equally illuminated field except directly opposite the 
extension of the luminous plane. Here if absorption is appreciable a vertical 
shadow will occur on the paper. No such shadow was obtained however. This 
supports the view that the stationary concentration of NH, radicals is very small 
and that the NH, decompose very quickly. The strength of the emission is 
not exceptional in view of the fact that the temperature at that point is about 
2 100%: 

Several additional facts can be gleaned from the analysis. For example, the 
strongest visible radiation comes from a point some distance from the position of 
maximum temperature where the reaction rate is highest. Also on the ammonia 
side the OH in absorption ceases at about 2,100°c., whereas on the O, side it 
is still very strong at that temperature. ‘This may be due to the lower partial 
pressure of OH in water vapour than in the HJO—O, mixtures, but it is rather 
striking that the OH ceases just where the NH, dehydrogenation begins. It may 
be an important dehydrogenation process, though it is unlikely to be the only one, 
otherwise the temperature distribution in the flame should be somewhat different. 
‘This may be shown by assuming that the OH does all the dehydrogenation since 
it is more exothermic in reaction with NH3, NH, and NH than the H or O atoms. 


Thus NH,;+OH —+NH,+H,0+~24 kcal. 
NH,+OH —+NH +H,0+-~24 kcal. 
NH +OH —+H,0+N_ +-.33 kcal. 

and N +N+M—-+N, +M* +225 kcal.* 


* The symbols M and M* are the usual convention for denoting collision with a third body 
tto take away energy of recombination N-+-N. 


PROC. PHYS. SOC. VOL. 62, PT. II—A (H. G. WOLFHARD AND W. G. PARKER) 


| 
OH 3064 a. OH 3428 a. 
Figure 4. 


800° c. 


600° c. 


BOY. 


200° c. 


DVaCs 


| | 
2200 2300 2400 


Absorption of NH; 
Figure 5. 


NH, 2245/2239 a. 


Figure 6. 


| | 
OH(00) NH OH(01) 


3064 a. 3360a. 3428 a. 
Figure 8. 


To face page 728 


Spectroscopic Examination of Flames at Normal Pressures 729) 


Hence a very big release of energy would occur at the point of dehydrogenation 
which does not accord with the temperature plot of the flame. It seems more 
likely that the product of dehydrogenation is mostly H, and not H,O and that 
the H, is subsequently burnt with oxygen and provides the greatest heat release 
of the reaction at about 5-6 mm. (Figure 7). 

A fuller investigation could be made with this technique and more quanti- 
tative data obtained from the band strengths would greatly help in completely 
elucidating the reaction mechanisms. 


§5. HYDROCARBON FLAMES 


Although no extensive work has yet been done with hydrocarbon flames, 
a few preliminary tests with methane-oxygen flames in this burner are worth 
recording. Onthe CH, side of the burner there was a narrow plane of strong 
luminosity due to radiant soot particles and on the O, side a weak blue continuum 
which was probably due to recombination of electrons and ions. These two 
vertical planes, visible to the naked eye, were apparently separated by a thin dark 
space. Figure 9 gives some of the results graphically represented as before. 


is ne 


° 
g000 C 


2000 °c 


> —4 
RG i 
N 


< OH ABSORPTION 
<e ‘ 
. | 


\ | 
INTENSITY 


\ OF BANDS 
% | 
\ (QuALitative ) 


CONTINUOUS 
EMISSION (SOOT ) 


19900 


\ Op 
“\ABSORPTION », 


ie) rd 4 G 8 16 
wy — 


EXTENSION THROUGH REACTION ZONE 


Figure 9. CH,-O, flat diffusion flame. 


The CH, diffuses from the left and the O, from the right; the distances in the 
flame are marked in millimetres on the abscissa. The visible continuous radiation 
(from carbon particles) extends from 2-8 to 4-0 mm. and the OH in absorption 
and emission just reaches the carbon region. ‘The OH bands are again much 
stronger on the O, side than at the corresponding temperatures on the CH, side. 
The maximum temperature (by Na line-reversal method) is close to the 
theoretical. It is also interesting to note that O, does not reach the luminous 
zone and the carbon particles cannot therefore burn in oxygen but probably react 
with OH or H,O. The region of greatest heat release is at 6:0 mm. and it seems 
very probable that this is due to combustion of H,—-O, and perhaps CO-O,. 

The information which this preliminary work gives relevant to carbon forma- 
tion may beindicated ; a discussion of this complex phenomenon will be presented 
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elsewhere. It is commonly supposed that carbon particles in the flame are 

formed by condensation of C, molecules or even C atoms which are sometimes 

observed in emission from premixed flames. This theory is quite untenable 

for the type of flame examined here because the C, is not observed prior to the 

_zone of carbon particles, i.e. on the fuel side of the luminous zone (Figure 9). 

Rather weak C, bands occur together with CH bands on the oxygen side of the 
zone. 

§6. CONCLUDING REMARKS 


It is hoped that sufficient evidence has been presented to show the potential 
value of the flat flame technique in the application of spectroscopy to diffusion 
‘flames. Very full information can be obtained in this way with fuels which 
-absorb light. For the study of hydrocarbon flames other than aromatic 
flames it may be necessary to use the near vacuum ultra-violet region, but this 
does not present very great difficulties. _Nitrogenous fuels and many other 
combustible inorganic compounds suchas compounds of sulphur can be very fully 
examined. The behaviour of fuel additives can be investigated and so also can 
‘the action of any promoters or inhibitors which give absorption spectra. The 
‘problem of carbon formation has already been mentioned and the technique of 
the flat diffusion flame is at present being used in this connection by the authors. 
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ABSTRACT. It is well known that some luminescent materials, notably zinc sulphide, 
‘show marked changes in dielectric constant and dielectric loss when excited. Recent work 
by the authors has established that the changes are due to the filling of electron traps, which 
have large effective diameter and high polarizability. 

The present paper describes measurements of the dielectric changes in phosphors 
‘containing two or more luminescence activators. Assuming that the dielectric changes are 
due to electron trapping, it is possible to derive information regarding the nature of electron 
traps and in particular their apparent association with the luminescence centres. The 
results support the contention that the introduction of the activating impurity into phosphors 
‘produces trapping states for excited electrons and that each trap forms part of a complex 
which includes the luminescence centre. 


SIN TROD LON 


ECENT studies of the extent to which retrapping occurs ia conventional 
phosphors show that it is a negligible process when the electrons are ejected 
from the traps by thermal energy only (Garlick aad Gibson 1948). This 

result cannot be reconciled with the conventional model of a luminescent material 
used by Randall and Wilkins (1945) and other authors. It was suggested, 
therefore, that each luminescence centre has one electron trap associated with it, 
‘the two together forming a closed system from which the excited electron did not 
normally escape. This view was supported by the experiments described by 
Garlick (1948) and by Garlick and Gibson (1949). Garlick showed that the 
introduction of impurities modified the trap distribution in the phosphor specimen. 
‘This was illustrated by the production of new peaks, due to the introduction of 
specific impurities, in the thermoluminescence curve. Garlick also showed that 
‘the thermoluminescence curves of phosphors containing two activators could be 
resolved into ‘colour peaks’, that is, the colour of the emission at a peak of the 
thermal curve glow depended upon the particular group of traps being emptied. 
‘The colour changes during thermoluminescence might be explained by migration 
of positive holes and electrons between different types of luminescence centres, 
as postulated in a recent theory due to Klasens et al. (Klasens, Ramsden and Chow 
‘Quantie 1948, Klasens and Wise 1948). However, this explanation is unlikely 
in many cases as the following experiment shows. If a phosphor with more than 
one type of luminescence centre is excited by radiation of suitable wavelength 
then it is possible to excite one type ot centre only. The resulting thermo- 
luminescence curve obtained after such excitation, using suitable optical filters, 
may be the same as that after excitation affecting both types of centre, provided 
that the same optical filters are used to select emission from the same luminescence 
centres. If this is so, then it is unlikely that any interchange of electrons or 
* Now at T.R.E., Malvern. 
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positive holes occurs between the different centres. In effect, sucha result supports 
the hypothesis mentioned above (Garlick and Gibson 1948) that retrapping of 
electrons does not occur in conventional phosphors, such as the zinc sulphides. 

Work on the increase in dielectric constant and loss which occurs in zinc 
sulphide on excitation (Garlick and Gibson 1947) has shown that the increase is 
due to the filling of electron traps. The present studies, described below, contain 
further evidence from dielectric measurements to show that electron traps and 
luminescence centres are closely associated. In particular, the effects of field 
frequency, temperature, exciting intensity and time of excitation have been 
studied. It has been shown that the dielectric constant change per filled trap is 
exponentially dependent on temperature, the exponent varying from phosphor 
to phosphor. The dependence on field frequency is typical of a dipole system 
having a relaxation time of about 10-7second. ‘The latter figure also varies 
among phosphor specimens of different constitution. 

It was noticed during the work on dielectric effects that both the relaxation 
time and the temperature exponent of the polarizability were, to a first approxi- 
mation, determined by the nature of the activating impurity. The latticestructure 
(cubic or hexagonal), the constitution of the bulk material (i.e. ZnS or mixed 
ZnS—CdS) and other parameters were of secondary importance. This suggested 
that the activating impurity producing the traps largely determined the dielectric 
properties of the traps. 


§2. EXPERIMENTAL RESULTS 
Garlick (1948) used the phosphor CaWO,-U as an example of a doubly 


activated material and gave the thermoluminescence curves of the material, using 
optical filters to select emission from the different types of luminescence centre. 

The same phosphor will be used as an example of the characteristics of the 
dielectric changes occurring when phosphors containing more than one activator 
are excited. 


(1) The Variation of the Dielectric Changes with Applied Field Frequency 


Measurements have been made by the methods described previously (Garlick 
and Gibson 1947) and over the same frequency range (75kc/s. to 20 Mc/s.). 
In a singly activated phosphor the relaxation time 7 is single valued or nearly so, 
anomalous dispersion of simple form occurring at about 2Mc/s. In a doubly 
activated phosphor, however, it is clear from inspection of the dispersion curve 
that there are two regions of anomalous dispersion and hence two values of the 
relaxation time. If the values are sufficiently different it is quite simple to treat 
the two dispersion regions separately and apply to each region the formula 
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given as equation (8) by Garlick and Gibson (1947)*. In Figure 1 the factor on 
the left-hand side of equation (1) is plotted against the square of the field frequency 
for the two dispersion regions considered separately. The slopes of the two lines 


* Ae is the change in the real part of the complex dielectric constant at the given field frequency, 
while Je) and Ae, are its values at zero and optical frequencies respectively. 
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give the two values of the relaxation time. The values obtained are 7, =31 x 10-8 
second, 7,=5-3x10-8second. Measurement of the variation in loss factor 
change Ae”, given by (Garlick and Gibson 1947) 
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shows that this factor has two frequency maxima corresponding to the two 
relaxation times already noted. As Ac” is a maximum when w7=1, values of + 
may be obtained. For the specimen of uranium activated calcium tungstate 
used they are 7,=23 x 10-8second and 7,=4:0 x 10-8second. The agreement 
with those deduced from the dispersion curve is as good as can be expected. 

In view of the results discussed above it seems reasonable to suppose that two 
different types of electron trap exist in this and similar materials. Each type of 
trap can be associated with one or other of the two types of luminescent centre in 
the material. In the above case of CaWO,-U we may talk loosely of ‘WO, 
traps’ and‘ Utraps’. It will be clear from the thermoluminescence curve given 
by Garlick (1948) that ‘U traps’ are deeper than ‘WO, traps’. 
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Figure 1. Variation of dielectric constant change with applied field frequency for a CaWO,-U 
phosphor. Curves plotted on assumption that phosphor has two types of polarizable centres, 
as explained in text. A: ‘WO, traps’; B: ‘ U traps’. 


4 
It is now important to ascribe the two values of 7 obtained by experiment to 


the ‘WO,’ and‘ Utraps’ respectively. One method is to determine the reciprocal 
of 7(=s, Garlick and Gibson 1947) for the two types of trap by luminescence 
measurements. ‘This can be done by the method due to Randall and Wilkins 
(1945), but high accuracy cannot be expected. Experiment indicates that the 
large value of 7 is associated with the deeper‘ Utraps’. This result was confirmed 
by a second method. If the decay of the dielectric constant change after 
cessation of excitation is measured at frequencies of the order of 5 Mc/s., the decav 
rate is fast, but becomes slower as the applied field frequency is reduced to a low 
value (~500ke/s.). At high frequencies the dipoles of small 7 are operative, and 
hence the smaller value of + must be associated with the shallower (in this case 
WO,,) traps. 
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In general, therefore, traps associated with given centres have a relaxation time 
+ peculiar to the impurity involved. Some typical results for zinc sulphide 
phosphors are summarized in the following table. 


Phosphor 7 (seconds x 10°) 


Zn traps Cu traps Mn traps 
ZnS-Zn 35 -- _- 
ZnS-Zn-Cu 34 i? — 
ZnS-Zn-Cu-Mn (1) 46 19 IDs, 
ZnS-Zn—Cu-Mn (2) 36 15 2c, 


(ii) The Variation of the Dielectric Changes during Thermoluminescence 


Measurements of the changes in dielectric constant and loss which occur 
during the emission of thermoluminescence have been made by methods already 
described (Garlick and Gibson 1947). The results obtained for doubly activated 
phosphors may be treated in the way described in the same paper. Ignoring the 
slow variation in polarizability with trap depth, the magnitude of the dielectric 
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Figure 2. ‘Temperature function f,.(T) for the dielectric constant change in a CaWO,-U phosphor 
derived from measurements of dielectric changes during thermoluminescence. 


constant change at any temperature T divided by the total number of filled traps 
remaining at that temperature is taken as a measure of the polarizability of the 
filled traps. Denoting the polarizability as f,.(7), it was shown by Garlick and 
Gibson (1947) that 
Tec Pye Const expe 1) oe) en ae ee (2) 

where « is a constant. 

The variation of f(T) with temperature for a doubly activated material 
(CaWO,-U) is shown in Figure 2, f,.(T) being plotted logarithmically. The 
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simple linear relation expected is not obtained, a marked inflection occurring at 
about 280°K. At about this temperature the emission changes from blue (WO, 
emission centres) to green (U emission centres). Apparently equation (2) is 
obeyed while ‘WO, traps’ predominate and while ‘U traps’ predominate, the 
inflection occurring intermediately. 

This result is explicable if it is assumed that the polarizabilities of ‘WO,’ and 
“U traps’ are different. As there is a drop in the ee ee 
curve at the inflection, it is clear that the polarizability of ‘WO, traps’ (Pywo,) 
is greater than the polarizability of ‘U traps’ (Py). The ratio P, 79,/Py can be 
found by extrapolating both curves to 0°K. and is 


Pwo,/Py =3-4. eevee (3) 


It will be clear that such a difference in polarizability will influence the decay 
of the dielectric constant change at room temperature, making the initial portion 
relatively faster. This is, in fact, observed. 

A difference in apparent polarizabilities of ‘WO,’ and ‘U traps’ is bound 
to arise, due to the difference in 7 value for each kind of trap. As the value of the 
constant « obtained from Figure 2 is almost the same for each kind of trap (indicated 
by the fact that the two linear parts of the curve ia Figure 2 are nearly parallel), 
it is likely that the difference in the values of « determines the difference in 
polarizability. The ratio of the apparent polarizabilities may be calculated from 
equation (1) as w, Ae, Ae,,, andzare known. It is found that, at the measurement 
frequency of 2:7 Mc/s., 

PwonlGeeo: | | | et ue vie (4) 


Agreement with the other experimental value above is better than might be 
expected. 


§3. DISCUSSION AND CONCLUSION 


Measurements of the increases in dielectric constant and loss of simple 
phosphors when excited suggest that the increases are due to the filling of electron 
traps (Garlick and Gibson 1947). If this result is assumed, dielectric measure- 
ments may be used to obtain more information on the nature of electron traps. 

The study of the dielectric changes occurring in doubly or triply activated 
materials shows that the theory already developed (Garlick and Gibson 1947) 
can only be extended to cover such materials if fuller assumptions are made. 
It is necessary to assume that the introduction of the impurity centres gives rise to 
the formation of electron traps, and further, that some measurable parameters 
of the filled traps are determined by the nature of the impurity. ‘This hypothesis 
has been suggested in previous work (Garlick and Gibson 1948, Garlick 1948), 
the measurable parameter concerned there being the depth of the electron traps. 
In the present case the parameter measured was relaxation time 7 of the filled 
electron traps. 

The relaxation time of such a dipole system is related to the diameter and other 
fundamental properties of the electron trap. The relation between 7 and the 
nature of the activating impurity thus gives further evidence in favour of the view 
that the impurity ion forms the electron trap, possibly by distortion of the local 
lattice. 
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Eigenfunctions following from Sommerfeld’s Polynomial 
Method 


By A. RUBINOWICZ 


Institute of Theoretical Mechanics, University of Warsaw 
MS. received 7th Fune 1949 


ABSTRACT. The eigenfunctions following from Sommerfeld’s polynomial method were 
determined in the general case. They are products of a convergence factor EF and a poly- 
nomial P given by a hypergeometric series. The form of E depends upon whether P is 
expressible by the ordinary or the confluent hypergeometric series. 


OMMERFELD’S polynomial method (Sommerfeld 1939) for the solution of 
the eigenvalue problems of the quantum theory is applicable only to a 


el 


» Se Ae ee ee 


limited group of such problems. It supposes that the eigenfunctions f(x) — 


are expressible in the form of a product of two functions f(x) = E(x) P(x). 
P(x) is a polynomial fulfilling a linear differential equation of the second order 
of the form 


7d a 
x?( A, + Box") ie ob 2x(A,+B,x") — (4, eBoe) P= Oe (1) 


where A,, B; are constants (A, 40) and h is a positive integer. 


The factor E(x) however is responsible for the convergence of the normalizatiom — 


integral and the fulfilment of boundary conditions. In a paper published. 
elsewhere (Rubinowicz 1949) it is proved that E(x) can appear only in two forms,. 
corresponding to B,0 and B,=0, provided that A,~0. If the differential 
equation of the eigenvalue problem has the form 


aC z) —4f +rApf=0, 
then we have 
(I) in the case B, 0, 
B= p-hy4l4a( Ao Bae) Pa R= ald ee tees (2 ay 
(II) in the case B, =0, 
JE =p txAiesp (17h) By A, ee (2b) 
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Since the differential equation of the polynomial (1) is reducible to the hyper- 
geometric differential equation (cf. Riemann-Weber 1901, p. 11) we can now write 
down the eigenfunctions for all cases accessible to Sommerfeld’s polynomial 
method. For this purpose we must transform equation (1) into the hyper-- 
geometric differential equation. Writing 


es ee | Be, (3) 
and 
PES ASN) Ba el a ce (4) 
we obtain the differential equation 
2 ae ese OP zs 
ee pag Ol aa Gye 8. enn (5) 
Here 
P2=Mg,  Py=M+2yMy, Pyo=Ny tpn +p(u—1)ng, (a 
pees a 
Q2=Mx, Gy =M,+2uIM, qgo=My tym, +u(u—1)me, | 
where 
WA, ty ti2A,+(h-1\A,t np =As, A 
my =h?By, m,=h{2B,+(h—1)Bz}, m=By | me 


Determining y in such a way that py vanishes, that is, from (6a) and (68), so 
that 
ph(wh —1)A, +2uhA, +A, =0, 


and introducing a new independent variable 


q2 By 
Z=-—-= ke ME Se, ie? be Suh ah (7 
, Pa” A, ) 
we get in place of (5) the hypergeometric differential equation 
a? d 
o(1—2) SE + fy—-(2 +B + 1)s} F - ap s=0. oe (8 a) 
‘The constants «, 8, y are here given according to (6 a) and (6b) by the equations 
ee eg ey Ce 
eee Seyi 52m 1) +241, freien (9a) 
anlar OT Sie eta MN ban ec 1) Bs 09 
ap = re pees p(w —1) PB, (uA(wh — 1) Bo + 2phB, + Bo}, 
ape ses (9d) 
Pi _ 4 l ( A, ) 9 5 
= == ee == — || PSs il +2u+1, yearns Ve 
fst ke A\” Ay \ 
From (9a) and (94) we obtain for €=« or =f the equation 
(uh — Eh)(wh — Eh —1)B, +2(uh—Eh)B,+By=0,  ...... (10) 


which is identical with Sommerfeld’s condition for the breaking off of the power 
series. 

The function ¢ which is a solution of the hypergeometric differential equation 
(8 a) is given by 
¢ = F(a, 8, y,2)= Fla, B,y, —(Be/Ae)x"], =... (11) 
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so that from (3), (4) and (11) we obtain for the polynomial P the expression 
P(x) = x** Fla, p, Y> =(5,/A.)0"1. AA dato (12a) 


In order that P(x) may be a polynomial, « (or 8) must be a negative integer. 

But we must remark that our transformation (7) can be used only in the case 
where B,-40, and is inapplicable if B, =0. 

In this second case we have to apply the substitution 


which transforms the differential equation (1) into 


pe es =i =e Se Sarees (8 5) 


The constant y is here determined by (9c), whereas « is given according to (10) 
by 


(8 b) is the differential equation of the confluent hypergeometric function F(a, y, 2). 
A solution of Sommerfeld’s differential equation (1) is therefore in case (II) given 
by 


Pye i (+. ae w) Je Poe Rea (12 d) 
1 


We are now able to write down the complete eigenfunctions f(x) = E(x) P(x). 

We get 
(1) in the case B,40 from equations (2), (9a) and (12a), 
fix) =p Bevenge—DA(] 4 Eater +D?2 Fig, By, —€) where &=(B,/A,)x"; 
(IL) in the case B, =0 from equations (26), (9c) and (12.4), 
B 
I(x) = prt Chaney Vie eci2 i (g) yw =) where C= : re xh, 
2 
The constants « or 8 are here negative integers. 

Using these expressions we can determine immediately the eigenfunctions in 
all cases accessible to Sommerfeld’s polynomial method. For this purpose we 
must know the values of A,, B; and h, which enable us to calculate the constants 
a, 8, y with the aid of the formulae given inthis paper. For all special cases we have 
dealt with in our former paper (Rubinowicz 1949) the values of A,, B; and h can 
be found there. 
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Origin of Low-Energy Electron Lines in the 
B-spectrum of **Au 


The low-energy part of the B-spectrum of 1°8Au was investigated, using the 254°-radial 
electrostatic field spectroscope that has been constructed in this laboratory. Several 
samples made from pure precipitated gold and also from gold foil were irradiated in the 
Harwell pile. It was found in all cases that two fairly intense conversion electron lines 
appear in the spectrum, one at energy 44 kev. and the other at 58 kev. Interpreted as due 
to internal conversion in gold in K and L shells respectively, they would correspond to y-ray 
energies of 0-125 Mev. and 0-072 Mev. The intensity of these lines was observed over 
five days and the half-life was found to be of the order of three days. 

The conversion-electron peak at 58 kev., which is broad, and due to 0:072 Mev. radiation, 
is in agreement with the observation of Jnanananda (1946), but there is a lack of agreement 
among observers about y-rays from 1°8Au, other than the well-known one at 0-4 Mev. 
Electron lines due to y-rays of very nearly the same energies as above have been reported by 
Valley (1941) and Helmholtz (1942) for !°7Hg, which decays by K-capture to *!®7Au, the 
y-rays, which are highly converted, being due to transitions from the excited states of 1®’Au 
to the ground state. It is suggested, therefore, that these conversion lines and the corre- 
sponding y-rays are to be associated not with the slow-neutron induced activity, !°8Au, but 
with the activity induced by fast neutrons during irradiation, giving !°’Pt by (n, p) reaction. 
eins, Aan (2) Peo i Aa. Excited, *°*Au produced -by the fadecay of "Pt 
would, on transition to the ground Jevel, emit the same y-rays, highly converted in gold, as 
are found in the K-capture decay of 1°’Hg. 

Krishnan and Nahum (1941), who produced 1°’Pt by (d, p) reaction on platinum, found 
that the 2-8 day isotope, as studied by absorption methods, emits strong discrete groups of 
low-energy electrons of secondary origin, the maximum energy being 0-120 Mev. Assuming 
K, L and M internal conversion in gold to be taking place, this would coirespond to a y-ray 
energy of 0-124 Mev., in good agreement with the energy of the y ray reported above. They 
also found evidence for a softer y ray from 1°’Pt and give the energy as 0-075 Mev., which 
again is in close agreement with the other y-ray energy mentioned above. The half-life for 
this isotope, 1°’Pt, as given by McMillan et a/. (1937) and Krishnan and Nahum (1941) is 
of the order of three days, which was also the value found for the decay of these conversion- 
electron lines. 

A long strip of gold foil, made into a small tight roll, was then irradiated at Harweli, so 
that the inner portion of the coil would be shielded from resonance neutrons, which produce 
the #*8Au activity. Measurement of the strength of 1°*Au activity as indicated by the 9-4 Mev. 
y-ray from equal areas of the shielded and unshielded portions of the foil showed thai the 
unshielded portion was substantially stronger than the shielded one, as is to be expected. 
On the other hand, the intensity of the convetsion-electron lines, relative to that of the con- 
tinuous f-ray background of 1°*Au, was greater in the spectrum of the shielded foil. This 
further strengthens the suggestion that the activity responsible for these y-rays and the 
conversion lines is induced in gold by fast neutrons and is not the '°%Au activity. 

A fuller account of this and other investigations of conversion lines from some other 
isotopes and also of the electrostatic B-ray spectroscope used for this purpose will be given 
at a later date. 

The author is very grateful to Professor P. B. Moon for his continued interest and 
guidance. He also wishes to thank Dr. D. M. Millest, with whom he collaborated in the 
design and construction of the f-ray spectroscope, and the Government of India for a 
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The Ultra-Violet Absorption Spectrum of SnO 


In the course of other work, we have recently obtained low-dispersion spectrograms of 
the absorption by SnO in the range 2040-3500 a. (see Plate). Two band-systems are known 
in this region: the less refrangible system, D—X, studied especially by Connelly (1933), 
and the E-X system, first obtained by Loomis and Watson (1934) in emission. The 
absorption plates add little to our knowledge of the D-—X system, but about twenty new 


The E-X System of SnO 
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‘bands have been assigned to the E-X system. Apart from the necessity of raising their 
v’-values by one unit, the vibrational analysis of Loomis and Watson is unchanged, and the 
‘main effect of our measurements, which are summarized in the Table, is to provide more 
information about the course of the upper-state vibrational levels. (The new values of 
AG’ are also more consistent than those of Loomis and Watson : the reason seems to be that 
for the vibrational analysis of a band-system of a molecule which gives an irresoluble 
complex of isotopic heads it is a positive advantage to use not too high dispersion.) 

In attempting to derive an expression for G,’, we used the values of G,” given by Jevons 
(1938) : G,”’=822-4(v’ +4) —3-73(v’ + 4)2.. No simple expression for G,’ was however found 
‘to fit the levels over the whole range, but up to v’=17 the equation 

G’=508-0(v’ + 4) —2:9(v’ + 3)?—5 x 10-9(v’ + 9)5 
is reasonably satisfactory. With this, ve=36,295 cm7!. 

The sharp drop in the vibrational interval beginning at v’=18 seems to be real and 
suggests that this point is not far from a dissociation limit. Graphical extrapolation gives 
Do’ ~ 1:2 ev., corresponding to a dissociation limit 5-6,ev. above z”=0. ‘This value of Do’ 
agrees fairly well with that suggested earlier by analogy with SnSe (Vago and Barrow 1946). 
Linear extrapolation of the ground-state vibrational intervals leads to Do’ =5-6 ev., so it is | 
possible that both states dissociate into the same atomic products. 

Band-systems of SnO are also expected at shorter wavelengths than the E—X system. 
‘On some of our plates weak red-degraded bands are to be seen with heads at about 2037, 
2050, 2054, 2069, 2073, 2084, 2102, 2117 and 2132 a. These may belong to a new 
‘system, but we have so far been unable to get pictures on which the bands are sufficiently 
well defined for accurate measurement. 

To complete the outline of the absorption spectrum of this molecule—so far as it is 
‘known at present—it may be mentioned that, contrary to statements in the literature 
(e.g. Sharma 1944), the stronger bands of the B and C systems (Pearse and Gaydon 1941) 
are also observed in absorption (Mrs. R. E. Richards, unpublished work): the vibrational 
-analysis of these bands is however still uncertain. 
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X-Ray Line Broadening in Metals 


It is the author’s belief that the diverse opinions of different workers on the explanation 
‘of x-ray line broadening in metals can be reconciled in terms of the hypothesis of Dehlinger 
and Kochendorfer (1939 a, b) and Kochendorfer (1944) that both particle size and stress 
effects are presented simultaneously. The theoretical work of Bragg (1942, 1949) gives 
‘mathematical support to the hypothesis. 

In the limiting case of pure particle size broadening the line breadth fp is related to the 
‘effective particle size ¢«, the x-ray wavelength A and the Bragg angle @ by the equation 
“Jones 1938) 

A 


ecos 0 


Bp= 


For pure stress broadening the breadth fg is related to the effective strain » (Stokes and 
Wilson 1944) and the Bragg angle by the equation 


Bg=7 tan 8, 


742 Letters to the Editor 


If both types of broadening are present, the resultant line breadth B should be obtained by- 
compounding the component lines as discussed by Jones (1938), Shull (1946) and Stokes. 
(1948), but it has been shown by Wood and Rachinger (1949) that for a commonly occurring- 
line shape the breadths are additive, ie. B=Bp+ Py. In this case we may write 


cos 6 1 sin 6 
sll BU Eo 


jae 
In a distorted lattice the effective particle size « must be interpreted as a measure of the- 
volume of regions in the lattice which diffract coherently. In anisotropic metals the strain 


Summary of line broadening measurements. 
. Aluminium filings, CuKa, Hall and Williamson (unpublished measurements). 
. Copper filings, ZnKa, Brindley (1939). 
. Copper filings, CuKa, Stokes, Pascoe and Lipson (1943). 
. Copper sheet, 20% reduction, CuKa, Dehlinger and Kochendorfer (1939). 
Copper sheet, 60%, reduction CuKa, Dehlinger and Kochendorfer (1939). 
. Copper sheet, 99:5% reduction, CuKa, Dehlinger and Kochendorfer (1939). 
. Rhodium filings, CuKa, Brindley (1938). 
. Rhodium filings, ZnKa, Brindley (1938). 
. Martensite rod, FeKa, Wheeler and Jaswon (1947). 


Key: 


distribution is not the same for different crystallographic planes, and a relationship of the- 
type 
Beos 6 4 20 sin 0 
4 = a 
A Gli 2X 
is frequently a better representation of the experimental breadths. In this equation Epjq- 
is the value of Young’s modulus for the direction perpendicular to the planes {hkl} and o- 


is the Laue breadth of the stress distribution function, which is assumed to be independent: 
of direction. 
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The accurate determination of the value of « from these equations is difficult, since the 
intercept which defines it depends largely on the breadths of lines for which sin 6/A is small. 
Metals rarely diffract at very low angles, and for such lines the total broadening is small and 
difficult to measure accurately. A number of published results have been analysed in this way 
and the best straight lines shown in the diagram determined by the method of least squares. 
It will be seen that the slopes of the lines (i.e. the stress effects) increase progressively from 
soft to hard metals. The intercepts are all positive with the exception of that for martensite, 
and a very small error in the steep slope of the graph would account for this physically unreal 
result. The mean value of all intercepts corresponds to «=10-' cm. Owing to the un-- 
certainties of the analysis no more than the order of magnitude of this quantity can be 
regarded as significant. This is about the same as that determined by Wood and Rachinger 
(1949) by a different method. The physical significance of this value may be understood 
from the theory of dislocations. There is abundant evidence (Taylor 1934, Brown 1941, 
Koehler 1941, Cottrell and Churchman 1949) suggesting that the maximum density of 
dislocations in most cold-worked metals corresponds to an average distance between 
dislocation lines of about 10-* cm. The coherently diffracting regions of the lattice ought to. 
have an effective particle size rather greater than this, since the dislocation lines are thought 
to be long compared with their separation. This suggestion agrees well with the experi- 
mental values, and it seems that direct calculations of line broadening based on the dislocation 
model should lead to a much clearer understanding of the problem. 
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Magnetic Viscosity in Mn-Zn Ferrite 


Ina previous paper (Street and Woolley 1949), referred to below as I, a theory of magnetic: 
viscosity, based on activation energy concepts, has been given. The relation between the 
rate of increase of the intensity of magnetization of a specimen and the time of observation 
equation (9) of I), is 

dl ipk 


af 
ares oe [exp (—Apgt)—exp (— Cd], 


where the symbols have the meanings given in I. The general solution of this expression is. 


: E 
AI=ipkT Eee Bit Chee a 2s Allee (1), 


\J being zero at time t=0. 
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In I it was found that only a limited range of t need be considered, viz. for values of t 
from t>te, where Ei(— Ct.) <Ei(—Apt,), to t<to, where Agty<1. Over this range (1) reduces 
to 

AI=ipkT log ¢-" constant. *7 "7 ete (2) 


‘The results of experiments made on magnetic viscosity in alnico satisfied this simplified — 
equation (2). However, in those cases where observations are extended over time periods | 
greater than the limits detailed above, it is necessary to use the exact form given in equation 
(1). A convenient test of the expression is provided by considering the experimental results 
for magnetic viscosity in Mn-Zn ferrite (Snoek 1947). The numerical values of the 
constanis of equation (1) were chosen to fit Snoek’s experimental results, and the equation — 
obtained was 


AB=0-875[Ei(—5-4 x 10-°t) —Ei(— 1-64 x 10-7) +5-72], .....- (3) 


where AB is the difference between the value of the magnetic induction of the specimen at 
any time t seconds and that at t=0. A plot of equation (3) and Snoek’s experimental 
results are shown in the Figure, and it will be seen that satisfactory agreement is obtained. 


41 AJ (gauss) 


1 10 100 1000S 10,000 100,000 
Time t (sec) 


Magnetic viscosity in Mn—Zn ferrite. 
The full line is calculated from equation (3). @ Experimental points obtained by Snoek (1947). 


The curve of best fit as drawn by Snoek is almost identical with the theoretical curve of 
equation (3). 
From the numerical equation, the following quantitative results may be obtained: 


E,/kT=5-72;  C=1-64x10- sec;  47ipkT=0-875. 


Thus, if it is assumed that J=290° k., then the upper value of the activation energy, Eo, is 
appropriate to a temperature of 1,660° k. and ip=1-72 x 10™ gauss/erg. This value of Ep is 
very much smaller than the activation energy suggested by Snoek, viz. 10,000° x. 

Equation (1) is of the same general form as that obtained by Becker and Dé6ring (1939) 
in their treatment of magnetic viscosity as the analogue of dielectric after-effect. However, 
the importance of equation (1) lies in the fact that the constants involved are expressed in 
terms of quantities which are characteristic of the domain processes which lead to magnetic 
viscosity. ‘Thus, although it is fairly certain that the domain processes are different for soft — 
iron, Mn—Zn ferrite and alnico, it is of interest to compare the available values of the va1ious 
constants of these three materials, assuming relation (1) to hold in each case : 


Res Mn-—Zn ferrite Alnico 
C (sec) 2>¢40"4 1:6x 10-7 == 
ip (gauss erg~?) 2:6 1038 171022 SSK 1 
Ey/k (° K.) 10,300 1,660 15,000 (estimate) 


* From analysis of results obtained by Richter (1937) using soft iron at 55° c. 
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Snoek consideis that the phenomena of magnetic viscosity and time deciease of perme- 
ability are two different aspects of identical domain processes. It follows therefore that the 
activation energy theory of magnetic viscosity given in I should account satisfactorily for 
time decrease of permeability. From preliminary experimental results using alnico 
specimens, it appears that the theoretical predictions are Soe Taree It is hoped to publish. 
further details of this work at a later date. 


The University, Nottingham. ; R. STREET. 
12th September 1949. J. C. WooLtLey. 


Becker, R., and Dorine, W., 1939, Ferromagnetismus (Berlin: Springer), p. 247 ff. 

RICHTER, G., 1937, Ann. Phys., Lpz., 29, 605. 

SNOEK, J. L., 1947, New Developments in Ferromagnetic Materials (Amsterdam: Elsevier), p. 55. 
STREET, R., and WooLLEy, J. C., 1949, Proc. Phys. Soc. A, 62, 562. 


REVIEWS OF BOOKS 


Scientific Foundations of Vacuum Technique, by Saut DusHman. Pp. xi+ 882. 
First Edition. (New York: John Wiley; London: Chapman and Hall, 
1949). 90s. 


The introduction of Gaede’s mercury diffusion pump marked the commencement of 
great advances in high-vacuum technique. A new industry—vacuum technology—has 
grown up, and a small-scale laboratory operation has been transformed within the past 
forty years into an industrial one, carried out on a scale that would have been almost 
inconceivable in the minds of the early workers in the field. Speed, and greater speed, 
is demanded for the exhaustion of all types of electronic equipment, cyclotrons, betatrons, etc.,. 
and for vacuum distillation, dehydration and the evaporation of metals and other substances. 
These higher speeds have been attained by, paying more attention to jet design in pumps. 

Only a few books have been written on the subject, and the first one to appear was that 
written by Dr. Dushman in 1922. The recent rapid advances have more than justified 
the present volume, which workers in this field have required for some time. In spite of 
its 882 pages, and of the fact that it contains a reference to the researches of practically every 
investigator in high-vacuum technique, it is no mere encyclopaedia, no catalogue of methods 
and results, but contains the theory which is so essential if successful operations are to be 
understood and appreciated. It is not only a complete revision of the earlier work, but it 
includes a very considerable amount of new material. The diagrams are numerous and 
good, and a useful feature is the copious references to original papers. 

The opening chapters deal with the kinetic theory of gases and their flow through tubes 
and orifices, features which play important parts in high-vacuum work and which con-- 
siderably affect pumping speeds. ‘The treatment of mechanical and vapour pumps is 
interesting and explicit. The réles-of diffusion and condensation are discussed, and the 
extreme importance of jet design is fully explained. Alexander’s modern type of pump is 
described and reasons are given why all vapour pumps do not function as diffusion ones— 
a conclusion which agrees with Langmuir’s view on their operation. Manometers receive: 
full treatment and the various indispensable leak detectors are included. ‘Then follows 
a description of other methods whereby gases may be removed from vessels, such as the 
sorption of gases and vapours by solids, and in particular by activated charcoal, silicates, 
cellulose and metals. Diffusion through the latter and degassing problems are also fully 
discussed. ‘The chapter on chemical and electrical clean-up of gases at low pressures gives. 
a comprehensive survey of the clean-up of gases by evaporated metals, technical getters,, 
incandescent filaments, cold and hot cathode discharges, and the production of extremely 
low pressures in sealed-off devices. The later chapters deal with vapour pressures, rates of 
evaporation and dissociation pressures—including useful information on vacuum distillation 
—and the deposition of films. 

The book is a complete treatment of the subject, one which should be in the hands of 
every worker in the fields of physics, engineering and chemistry. 

F. H. NEWMAN- 
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Thermocinetique, by PIERRE VERNOTTE. Pp. xxi+459. Publications Scientifiques 
et Techniques du Ministére de I’Air. (Paris: Service de Documentation et 
d’Information Technique de |’Aéronautique, 1949.) No price. 


3, 


Monsieur Vernotte is a well-known writer on thermal subjects, and those familiar with 


his work are aware that he does not blindly follow convention without asking the reasons 
for it. In this book, which deals with the subject that is sometimes described as Heat 
Transfer, he opens his challenge to tradition by objecting to the traditional names, and 
selecting eventually that of thermokinetics, so that, between them, thermokinetics and 
thermodynamics cover the main part of the subject of Heat. 

After this, and before we begin on the main subject, we find a stimulating discussion 
of the question whether the unit of heat is indeed a unit of energy, and whether, therefore, 
many of us are right in saying that the erg or joule is the proper unit for heat quantities; 
the author thinks it is not, though it may be used for convenience. 

The main part of the book deals with conduction, convection and radiation. In 
dealing with the first of these, Vernotte sets up the differential equation in the usual way, 
with a perfectly correct warning that it is only strictly applicable if the thermal constants 
do not vary with temperature. He claims also that the passage from a difference between 
two first derivatives to a second derivative is not really valid, though he accepts it in 
practice. The reviewer confesses to an inability to understand the grounds on which the 
objection is based. : 

Having set up the equation, he solves a number of problems; these are well chosen 
for their illustrative value or for practical use, but they do not illustrate as many methods 
as might perhaps have been expected in a book of this size. There is no reference to the 
value of conjugate function methods in two-dimensional problems, nor of the Laplace 
transform. 

Convection is interestingly treated, and the main laws of radiation are brought out. 

Altogether, as mentioned earlier, the value of the book is as a stimulus to our own 
thinking, rather than as a compendium of methods or results, but for its actual purpose 
it is an extremely successful effort. iesass 


CORRIGENDUM 
“The Physical Basis of Life”, by J. D. BERNAL (Proc. Phys. Soc. A, 1949, 62, 537). 


Page 540, last clause of footnote *, viz.: ‘‘ and secondly, they must originate and develop 
out of some pre-existing system, or, in plain English, they must work and they must have 
got there in the first place ’’ skould occur in the text immediately after the words ‘‘ dynamic 

stability ”’ (line 22). 
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ABSTRACTS FOR SECTION B 


An Automatic Recording Apparatus for the Study of Flow and Recovery in Metals, 
by E. N. pa C. ANDRADE and A. J. KENNEDY. 


ABSTRACT. An apparatus is described which records continuously on photographic 
paper the extension-against-time curve of a metal wire creeping under stress. The length 
and time scales are recorded on the paper at the same time as the extension, so that deforma- 
tion of the paper during development and drying does not affect the accuracy of the record, 
which allows an extension of up to 10 cm. to be read to within 0:02 mm. The apparatus also 
automatically removes and restores, repeatedly if desired, the load at times which can be set 
before the start of the experiment. 


An Investigation of the Mechanical Properties of Materials at very High Rates of 
Loading, by H. Ko.sxy. 


ABSTRACT. Amethod of determining the stress-strain relation of materials when stresses 
are applied for times of the order of 20 microseconds is described. "The apparatus employed 
was a modification of the Hopkinson pressure bar, and detonators were used to produce 
large transient stresses. ‘Thin specimens of rubbers, plastics and metals were investigated 
and the compressions produced were as high as 20 % with the softer materials. It was found 
that whilst Perspex recovered almost as soon as the stress was removed, rubbers and poly- 
thene showed delayed recovery, and copper and lead showed irrecoverable flow. The 
phenomenon of delayed recovery is discussed in terms of the theory of mechanical relaxation 
and memory effects in the material. 
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A Method of Measuring some Dynamic Elastic Constants and its Application to 
the Study of High Polymers, by K. W. HILuirr. 


ABSTRACT. The preliminary work on the measurement of velocity of propagation and. 
attenuation of longitudinal sound oscillations in high polymer filaments already described 
by Hillier and Kolsky has been continued. It was concluded in that paper that temperature 
control would be advantageous, and that by varying the temperature and the frequency of 
oscillations more useful information could be acquired. An apparatus is described whereby 
the velocity of the propagation and attenuation of longitudinal oscillations throughout the — 
range 0° to 50° c. and 500 c/s. to 30 kc/s. can be measured. The results obtained with one 
material, polythene, are discussed in terms of several theories of the elasticity of high 
polymers, and the constants of the several equations considered are calculated from the data. 


The Theory of the Microscope—IV: The Boundary-Wave Theory of Image 
Formation, by L. C. Martin. 


ABSTRACT. The theory of boundary waves (Young, Rubinowicz, and others) is developed ~ 
for the two-dimensional case and applied to discuss image formation for small apertures, 
obstacles and phase-retarding laminae. The theoretical predictions are illustrated by 
diffraction photographs. It is hoped thus to contribute towards a fullér understanding of 
microscope images for bright-field, dark-field and phase-contrast conditions. 


Wave-front Aberrations of Oblique Pencils in a Symmetrical Optical System = 
Refraction and Transfer Formulae, by W. WEINSTEIN. 


ABSTRACT. ‘The wave-front aberrations which occur in a pencil traversing a symmetrical 
optical system at a finite field angle are defined ina new way. Refraction and transfer 
formulae are derived for all the aberration coefficients up to and including those depending 
on the fourth power of the aperture of the pencil. 


The Temperature Variation of the Long-Wave Limit of Infra-Red Photo- 
conductivity in Lead Sulphide and Similar Substances, by 'T. S. Moss. 


ABSTRACT. Measurements have been carried out on the spectral distribution of the 
photoconductive effect in layers of lead sulphide, selenide and telluride, at temperatures 
ranging from room temperature to that of liquid hydrogen. It is shown that cooling has 
a marked effect on the long-wave limit of sensitivity, which increases by as much as two 
microns on cooling with liquid hydrogen. 

Brief descriptions of spectral measurements on these materials have already been published 
by the author. It is shown that depopulation of energy levels due to thermal activation of 
electrons is not responsible for the shift of spectral sensitivity. "The suggestion is made that 
the shift may be related to the variation of dielectric constant with temperature. 
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‘THE PROCEEDINGS OF THE PHYSICAL SOCIETY 


-HYSICAL SOCIETY SPECIALIST GROUPS 


PTICAL GROUP 


-'The Physical Society Optical Group exists to foster interest in and development of all branches of 
ptical science. ‘To this end, among other activities, it holds meetings about five times a year to 
iscuss subjects covering all aspects of the theory and practice of optics, according to the papers offered. 


OLOUR GROUP 


The Physical Society Colour Group exists to provide an opportunity for the very varied types of 
rorker engaged on colour problems to meet and to discuss the scientific and technical aspects of their 
rok. Five or six meetings for lectures and discussions are normally held each year, and reprints of 
apers are circulated to members when available. A certain amount of committee work is undertaken, 
nd reports on Defective Colour Vision (1946) and on Colour Terminology (1948) have already been 
ublished. 


OW TEMPERATURE GROUP 


The Low Temperature Group was formed to provide an opportunity for the various groups of 
eople concerned with low temperatures—physicists, chemists, engineers, etc.—to meet and become 
amiliar with each other’s problems. ‘The Group seeks to encourage investigations in the low temperature 
eld and to assist in the correlation and publication of data.’ 


COUSTICS GROUP 


The Acoustics Group was formed to meet the long felt need for a focus of acoustical studies im 
Sreat Britain. The scope includes the physiological, architectural, psychological, and musical 
ispects of acoustics as well as the fundamental physical studies on intensity, transmission and 
ibsorption of sound. The Group achieves its object by holding discussion meetings, by the circula- 
jon of reprints and by arranging symposia on selected acoustical topics. 


Further information may be obtained from the Offices of the Society ; 
1 LowTHER GARDENS, PRINCE CONSORT Roap, LONDON S.W. 7. 


BULLETIN ANALYTIQUE 


Publication of the Centre National de la Recherche Scientifique, France 


The Bulletin Analytique is an abstracting journal which appears monthly in two parts, Part I 
covering scientific and technical papers in the mathematical, chemical and physical sciences and their 
applications, Part II the biological sciences. 

The Bulletin, which started on a modest scale in 1940 with an average of 10,000 abstracts per 
part, now averages 35 to 40,000 abstracts per part. The abstracts summarize briefly papers in scientific 
and technical periodicals received in Paris from all over the world and cover the majority of the more 
important journals in the world scientific press. The scope of the Bulletin is constantly being enlarged 
to include a wider selection of periodicals. 

The Bulletin thus provides a valuable reference book both for the laboratory and for the individual 
research worker who wishes to keep in touch with advances in subjects bordering on his own. 

; A specially interesting feature of the Bulletin is the microfilm service. A microfilm is made of 
each article as it is abstracted and negative microfilm copies or prints from microfilm can be purchased 
from the editors. 

The subscription rates for Great Britain are 4,000 frs. (£5) per annum for each part. 
Subscriptions can also be taken out to individual sections of the Bulletin as follows : 


frs. 

Pure and Applied Mathematics—Mathematics—Mechanics 550 14/6 
Astronomy—Astrophysics—Geophysics .. 700 18/- 
General Physics Thermodynamics—Hleat-Optics—Blec- 

tricity and Magnetism .. , : 900 22/6 
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General Chemistry—Physical Chemistry .. 325 8/6 
Inorganic Chemistry — Organic res — “Applied 

Chemistry—Metallurgy : 1,800 45/- 
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Biochemistry—Biophysics—Pharmacology ai 900 22/6 
Microbiology—Virus and Phages .. ats ae 600 15/6 
Animal Biology—Genetics—Plant Biology _ ae oe 1,800 45/- 
Agriculture—Nutrition and the Food Industries . . we 550 14/6 


Subscriptions can be paid directly to the editors : Centre National de la Recherche Scientifique, 
18, rue Pierre-Curie, Paris 5tme. (Compte-chéque-postal 2,500—42, Paris), or through Messrs. H. K. 
Lewis & Co. Ltd., 136, Gower Street, London W.C. 1. 
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